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15977 SOIL-PIPE INTERACTION OIL PIPELINES 


KEY WORDS: Corrosion; Pipelines; Powerplants; Protective coatings; Soil 
mechanics; Soil-structure interaction 


ABSTRACT: study presented analyzing the distress frequently encountered 
buried heated oil pipelines due the expansion and contraction producing soil-pipe 
interface shear and cyclic movements. This movement can damage the coating the 
pipe and result corrosion failure. The magnitude the cyclic shear stress and 
movement depends temperatures, soil properties, and pipeline design items such 
the length casing, bend, and straight pipe sections. Two important factors are the 
coefficient horizontal subgrade reaction and the soil-pipe interface shear parameter. 
The interface shear parameter unique the heated pipeline problem and can serve 
coating performance and selection criteria. Methods calculate pipe movements 
and ways restrain such movements are analyzed. Actual field measurements 
heated oil pipeline are compared with theoretical calculations. 


REFERENCE: Yen, Bing C., Tsao, Ching H., and Hinkle, D., “Soil-Pipe 
Interaction Heated Oil Pipelines,” Transportation Engineering Journal, ASCE, Vol. 
107, No. Proc. Paper 15977, January, 1981, pp. 1-14 


15987 HIGH HEAD PERFORMANCE CONCRETE PIPE 


KEY WORDS: Bells; Bending stress; Concrete pipes; Design; Diameters; 
High head; Performance; Pressure; Prestressing; Ring theory 


ABSTRACT: This paper deals with the manufacture and subsequent testing 600 psi 
(4137 kPa) 189-in. (4,800-mm) inside diamater prestressed concrete pipe with four 
layers prestressing wire. Existing barrel design theory and newly developed theory 
for determining joint ring stresses are evaluated and shown correlate with test data. 


REFERENCE: Davidson, James W., “High Head Performance Concrete Pipe,” 
Transportation Engineering Journal, ASCE, Vol. 107, No. Proc. Paper 15987, 
January, 1981, pp. 15-23 


15993 PERMEABILITY LIME STABILIZED SOILS 
KEY WORDS: Clayey soils; Curing; Density; Lime; Permeability; Soil lime 


ABSTRACT: The permeability two soils stabilized with lime when compacted 
various molding water contents has been investigated the laboratory. poorly 
graded river sand and sandy silty clay were used. The effect lime content, curing 
age, dry density, and type soil the permeability were investigated. The addition 
lime reduced the permeability the investigated sand various moisture contents, 
while the permeability the clay increased when compacted the optimum moisture 
content. The permeability all mixtures was found decrease with curing age. There 
existed molding water content for the sand-lime mixture which the permeability 
reached minimum. The permeability the investigated sandy silty clay increased 
with increasing lime content the same dry density. 


REFERENCE: El-Rawi, Nagih M., and Awad, Amir A.A., “Permeability Lime 
Stabilized Soils,” Transportation Engineering Journal, ASCE, Vol. 107, No. TE1, Proc. 
Paper 15993, January, 1981, pp. 25-35 
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15994 TAX ALLOCATION STEPWISE REGRESSION 


KEY WORDS: Alligatorweed; Distribution patterns; Regression analysis; 
Revenue; Statistical tests; Taxes; Transportation 


ABSTRACT: regression analysis model has been developed derive alternate 
distribution formulas for statewide fuel tax revenue. The model combines predictor 
variables such road mileage, population and vehicle registrations way that best 
fits the current needs local jurisdictions. analysis using data for Alabama’s 
counties produced distribution formulas that exhibited significant improvement over 
the formula use the present time. examination statistical tests assess the 
adequacy alternate distribution formulas also presented. 


REFERENCE: Bell, Lansford C., and Shelton, Russell D., Allocation Formulas 
Stepwise Regression,” Transportation Engineering Journal, ASCE, Vol. 107, No. 
Proc. Paper 15994, January, 1981, pp. 37-44 


15982 SOIL-PIPE INTERACTION AND PIPELINE DESIGN 


KEY WORDS: Deformation; Earth pressure; Failure; Finite element 
method; Pipelines; Soil deformation; Steel pipes 


ABSTRACT: The problem steel pipe bend design with and without lateral soil 
restraint studied, and the effect backfill soil deformation the bend design 
postulated. Studies show that finite element model utilizing nonlinear soil properties 
superior conventional earth pressure theories and can result economic 
design. The effect backfill type and the depth soil cover resisting the outward 
radial forces bend locations are examined. Practical design and construction 
recommendations are presented. 


REFERENCE: Ahmed, Syed, McMickle, Ralph W., and Brassow, Carl L., “Soil-Pipe 
Interaction and Pipeline Design,” Transportation Engineering Journal, ASCE, Vol. 107, 
No. Proc. Paper 15982, January, 1981, pp. 45-58 


15968 SUPERPLASTIC CONE-SHAPED SANDWICH CORES 


KEY WORDS: Aerospace engineering; Cones; Cores; Elastic properties; 
Energy methods; Engineering mechanics; Projection; Sandwich structures; 
Shear modulus; Structural design 


ABSTRACT: new sandwich core formed using the superplastic sheets and the 
vacuum pressure- gas pressure-forming technique, both. capable providing 
opportunities for structural optimization and better bonding with the facing sheets than 
the conventional honeycomb core. The shear modulus truncated, hollow, cone- 
shaped core with square symmetry determined. Energy principles are used for 
finding the upper and lower bounds solution. Results presented graphic form 
indicate that these bounds are indeed very close. The methodology developed here 
effective for providing excellent engineering solution the problem. The new 
superplastically-formed core also offers competent approach for those who seek 
more efficient metallic structural particularly when used combination 
with the difussion bonding technique for weight and cost savings. 


REFERENCE: Ueng, Charles E.S., “Superplastic Cone-Shaped Sandwich Cores,” 
Transportation Engineering Journal, ASCE, Vol. 107, No. Proc. Paper 15968, 
January, 1981, pp. 59-68 
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15984 IMPACT BEHAVIOR BARRIERS NONLEVEL TERRAIN 


KEY WORDS: Automobiles; Barriers; Guardrails; Highways; Impact 
severity; Safety; Slopes; Terrain; Tests; Vehicle trajectories 


ABSTRACT: Six full-scale crash tests were conducted evaluate the impact behavior 
two widely-used roadside barriers when placed 6:1 side slope. Four the six 
tests involved standard W-beam rail metal post G4(1S) system< and the other 
two tests involved three-cable barrier mounted metal posts (G1 system). The tests 
were conducted and evaluated accordance with national guidelines for testing 
roadside appurtenances. Based the results, recommended that the G4(1S) 
system not placed 6:1 steeper side slopes. Placement the system 
6:1 side slope acceptable provided ample space exists between the barrier and the 
hazard shielding. 


REFERENCE: Ross, Hayes E., Jr., and Smith, Darrell G., “Impact Behavior 
Barriers Nonlevel Terrain,” Transportation Engineering Journal, ASCE, Vol. 107, 
No. Proc. Paper 15984, January, 1981, pp. 69-79 


15996 PROCEDURES FOR ESTIMATING TRIP ATTRACTIONS 


KEY WORDS: Forecasting; Land use zoning; Linear regression analysis; 
Models; Transportation planning; Trip generation; Urban transportation 


ABSTRACT: For many large metropolitan areas, gross land area, stratified broad 
land-use categories, such residential, commercial, industrial, and open space, 
usually the only readily obtainable data outside the central business district. Need 
exists for reliable techniques that can used for trip-attraction estimation for such 
urban areas. This paper presents set equations that have been calibrated, using 
multiple linear regression analysis, compute zonal person-trip attractions for various 
trip purposes function gross land area. For each trip-purpose, equations have 
been further stratified pseudo-employment density ranges and dummy variables that 
reflect the presence major trip attractors capable satisfying the trip purpose 
question. For zones falling within very high density areas, equations have been 
calibrated function zonal gross floor area. Equations have alsc been developed 
for estimating zonal nonhome-based and truck trip ends. 


REFERENCE: Kumar, Ashok, “Simplified Procedures for Trip 
Attractions,” Transportation Engineering Journal, ASCE, Vol. 107, No. Proc. 
Paper 15996, January, 1981, pp. 81-99 


U.S. Conversion Factors 


accordance with the October, 1970 action the ASCE Board Direction, which stated 
that all publications the Society should list all measurements both U.S. Customary and 
(International System) units, the following list contains conversion factors enable readers 
compute the unit values measurements. complete guide the system and its 
use has been published the American Society for Testing and Materials. Copies this 
publication (ASTM E-380) can purchased from ASCE price $3.00 each; orders must 
prepaid. 

All authors Journal papers are being asked prepare their papers this dual-unit format. 
provide preliminary assistance authors, the following list conversion factors and guides 


are recommended the ASCE Committee Metrication. 


convert 


inches (in.) 

feet (ft) 

yards (yd) 

miles (miles) 

square inches (sq in.) 
square feet (sq ft) 

square yards (sq yd) 

square miles (sq miles) 
acres (acre) 

cubic inches (cu in.) 

cubic feet (cu ft) 

cubic yards (cu yd) 

pounds (Ib) mass 

tons (ton) mass 

pound force (Ibf) 

kilogram force (kgf) 

pounds per square foot (psf) 
pounds per square inch (psi) 
U.S. gallons (gal) 

acre-feet (acre-ft) 


millimeters (mm) 
meters (m) 

meters (m) 

kilometers (km) 

square millimeters 
square meters 
square meters (m’) 
square kilometers (km’) 
hectares (ha) 

cubic millimeters 
cubic meters 

cubic meters 
kilograms (kg) 
kilograms (kg) 

newtons (N) 

newtons (N) 

pascals (Pa) 

kilopascals (kPa) 

liters (L) 

cubic meters 


Multiply 


25.4 
0.305 
0.914 
1.61 

645 
0.093 
0.836 
2.59 
0.405 

16,400 
0.028 
0.765 
0.453 

907 
4.45 
9.81 

47.9 
6.89 
3.79 

1,233 
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INTERACTION HEATED PIPELINES* 
(Reviewed the Pipeline Division) 


INTRODUCTION 


Many utility companies operate buried pipeline networks for transporting heated 
fuel oil from tank farms power generating plants. Depending the demand 
generating stations and requirements for storage capacity, the heated oil 
transported intermittently between storage tanks and generating stations about 
typical section pipeline system. The intermittent flow operation 
requires that the fuel oil preheated before transporting. The intermittent 
operation also causes periodic pipeline thermal expansion and contraction and 
produces soilpipe interface shear and cyclic movements which may result 
damaged pipelines. High cost maintenance efforts have been necessary recent 
years repair and replace damaged pipelines. 

Fig. shows typical failure pipe and its protective coating. The 
coating suggests that the thermal expansion and contraction induced excessive 
movements and cyclic shear stresses between the pipe and the surrounding 
soil. this case, the are orientated perpendicular the longitudinal 
axis the pipeline. Fig. the orientation wrinkles near the pipe bend 
suggests lateral and radial movements the pipe. Figs. and where 
the coatings are wrinkled and longer protect the pipe, corrosion-caused pipe 
pitting has started. Fig. shows pitted areas steel pipe that failed either 
service under hydrostatic test. 
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Although considerable research has been conducted the soil-pipe interaction 
problem buried pipelines, only limited nonproprietary literature available 
heated oil pipelines. state-of-the-art survey available and references, 
thereto), including survey existing computer programs (6). However, 
sophisticated computer analyses are often limited the imprecise and sometimes 
poorly understood soil parameters. Furthermore, often impractical apply 


Note wrinkle orientation 
change near pipe bend 


FIG. Pipe Coating Caused Excessive Movement Pipe Sur- 
rounding Soil 

the same level effort for the analysis used for pipelines nuclear power 
plant heated-soil pipeline which may span many kilometers length through 
various soil conditions. Therefore, simplified approach compatible with the 
unique nature the thermal-induced movements heated oil pipelines has 
been called for the design engineer (1) and coating manufacturers (4) 
recent years. 


= 


HEATED OIL PIPELINES 


this article, the writers will: (1) Characterize the heated oil pipeline problem 
from geotechnical engineering viewpoint; (2) analyze soil parameters pertinent 


Note wrinkles in radial 
and lateral directions 


FIG. Pipe Coating Radial Direction Pipe Bend (note pitting 
exposed pipe) 


4» 


FIG. 3.—Pipe Failed Casing Entrance (note excessive pitting exposed pipe) 


simplified soil-pipe interaction model for analysis and design; and (3) consider 
examples and compare field measurements with the theoretical analysis. 
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From geotechnical engineering viewpoint, heated oil pipelines are charac- 
terized the following special features: 


Potentially, large force can exerted soils pipe bends. Due 
the relatively large temperature differential between intermittent operations 


the pipeline and the high internal operating pressure, the pipeline has potential 
exert large forces the surrounding soil pipe bends, e.g., assuming the 


-T2=4T 


FIG. 4.—Notations Pipe Bend 


Pian View 


FIG. 5.—Soil Bearing Area Pipe Bend 


backfill bends densely compacted that the pipe can considered 
fully restrained (Fig. 4), the maximum axial force resisted the 
soil bends may estimated by: 


max 


which the maximum force resisted the soil the bend 
assuming nonyield soil; area the steel pipe section; modulus 
elasticity pipe; thermal expansion coefficient pipe; temperature 
differential, degrees; Poisson’s ratio pipe material; average hoop 
stress; and pipe wall thickness. For instance, for typical 
16-in. (406 mm) steel pipeline with wall thickness in. mm) operated 


I cn N. aa 
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temperature differential 140° (60° and internal pressure 1,000 
psi (6.9 the maximum axial force 361 kips (1.6 MN) 45° 
bend. This large force resisted the soils pipe bends. 

Note also that Eq. independent pipe length, and the maximum force 
exerted the soil depends only operating temperature, pressure, and the 
geometry the bend. Actually, the forces transmitted the soil are lower 
because the soil the bend not infinitely rigid and does not fully restrain 
the pipeline; the actual force depends the resistance soil deformation 
pipe bends. The stiffer the soil the bend, the larger the force will be. 
But for the time being, suffice say that the heated oil pipeline capable 
exerting large forces the soil bends and these forces increase significantly 
the operating temperature increases. Under such large forces, the soil 
likely yield causing movements between the pipe and soil. This leads 
the second special feature. 

Cyclic Movements.—The cyclic movements induced intermittent opera- 
tions occur locations where the pipeline changes direction (bends), changes 
size, terminates. These movements also cause shear stresses soil-pipe 
interface. these movements are not sufficiently restrained and the shearing 
stresses become excessive, the coating may damaged. Methods restrain 
the pipe movements are commonly the use large pipe bend radius installation 
large thrust blocks the bend. These are expensive methods and frequently 
use large pipe bends infeasible will examined next. 

Limitation Use Large Pipe Bend Radius.—Using thrust blocks 
large radius bends restrain the pipeline movements may subject practical 


limitations. For pipe bends, the effective soil bearing area, may estimated 
(Fig. 5): 


which ratio bend radius pipe radius; and all other notations have 
been defined before. 

given street corner right-of-way, the pipeline designer often limited 
how large bend can used. small value means small soil area 
restrain pipeline movement. Therefore, from designer’s viewpoint 
essential provide quantitative technique estimate the pipe movements 
terms and soil properties. 

Cased Pipe Crossing.—Casing are generally used locations where pipelines 
cross under roadways. Pipelines are supported low-friction spacers within 
the casing, and therefore the pipe not restrained the surrounding soil. 
Movement and potential coating damage occur outside the casing where the 
surrounding soil contact with the pipe, e.g., Fig. shows pipeline that 
failed near casing entrance. Therefore, essential for the engineer 
able estimate the amount pipe movement the casing entrance. 


Although buried pipelines may have configurations consisting bends, casings, 
straight sections, and combinations thereof, these configurations can de- 
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composed into two generic elements shown Fig. combination these 
generic elements can simulate most pipeline configurations shown Fig. 
Thus, other pipeline design configurations can calculated once the mathe- 
matical model the two generic elements are solved. 

the mathematical models the pipe bends have been replaced soil 
restraining force For the straight pipe section where the soil contact 
with the surrounding soil, there soil-pipe interface shearing force acting 
against the pipe. Within the casing, there contact between the pipe and 


Section zero pipe movement 


FIG. 6.—Physical and Mathematical Models: (a) Generic Element and (b) Generic 
Element 


Analysis Generic Element Straight Pipe and Bend.— When the pipe heats 
up, the pipe expands. The soil-coating interface shear force acts along the 


pipe resist expansion and the pipe bend counters the expansion with force 
The stresses the pipe wall 


qdx 


(a) 
Section zero pipe movement 
q 
(b) 


TE1 HEATED OIL PIPELINES 
distance from the section zero pipe movement 


and the pipe strain the summation thermal-induced strain plus the strain 
caused soil stresses, i.e.: 


=— = — + (a: -v-o 


which pipe movement section Let C-u, i.e., assuming the 
soil-pipe interface shear force proportional the pipe movement. can 


FIELD CASE 


0 
Element A connect 
to element A. 


Element B connect 


to element B. 


Element B connect to 
element B and x = = 


(see Figure 6). 


Element B connect 
to element A. 


FIG. Showing Pipeline Configurations Combining Generic Elements 


shown that the differential equation describing the relationship pipe 
movements 


which the soil-pipe interaction parameter relating pipe property the 
soil-pipe interface shear: 


? 
@ 
0 
7 
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The physical significance and will considered later. The general 


Bx —Bx 


Parameters and may solved using boundary conditions: (1) 
which coefficient horizontal subgrade reaction; and pipe movement 
can shown that the pipe displacement given by: 


The stress the pipewall given 


1 
(10) 


Ss 


The interface shear given 


Analysis General Element Casing and Bend.—The derivations differential 
equations for generic element [Fig. identical generic element 
Eqs. 4-8, except for boundary conditions. can shown that 


The two unknowns, and C,, can obtained solving Eq. and Eq. 
simultaneously. Once and are known, the pipe movement 


TE1 
Cc 1/2 
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The soil-pipe interface shear distance from the casing entrance is: 

the pipe axial stress, distance from the casing entrance is: 


representing the relative stiffness the soil-coating-pipe system. large 
value means strong soil-coating-pipe bond that capable providing large 
restraining force along the pipe and limiting the pipe movement. small 
value represents weak soil-coating-pipe bond which provides little restraining 
force limit pipeline movements. show the significance B1, Fig. shows 


x 
0.4 


ad 


Pipe Length Ratio, x/! 


FIG. Soil-Pipe Property Pipe Movements for Pipeline with 


simulates condition such where pipeline terminates with 
bend, the soil the bend very soft loose such that offers 
resistance, both, could condition such that there gap between 
pipe bend and soil left from the previous cycle intermittent pipeline operation. 
can seen that for strong soil-coating-pipe system, pipe movements will 
small but most the movements occur within 20% the pipe length near 
the pipe end. But, because the high stiffness, the shear force within this 
20% length zone can too great resisted the adhesion between pipe 
and coating, too large resisted the shearing strength the coating. 
Pipe damage may ensue. the other hand, for soft soil-coating-pipe system 
low values, the shearing force along the pipe will small. But, because 
the low shearing force, the pipe movements will large and more evenly 
distributed along the pipe shown Fig. Depending the property 


. at 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 . 1.0 
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the coating used, the movements may too large for the coating resist 
and damage may ensue. Therefore, (and thus the value) unique 
criterion measure the performance and behavior the soil-coating-pipe 
interaction hot oil pipelines. 

Location Section Zero-Pipe-Movement.—Section 0-0 Fig. represents 
location which the pipe has movement relative the soil, 
Physically, Section 0-0 could the location: (1) Where the pipe connected 
rigid structure; (2) symmetrical condition with the existence identical 
bend configuration the left Section 0-0. Section 0-0 may located for 
straight pipe connecting different bends which pipe movements each 
bend are known: and pipeline length between 
the two bends. the soil conditions and bend geometry are different, 
iterations will required locate Section 0-0 using initially assumed displace- 
ments each bend. 

Influence Pipe Bend Radii and Coefficient Soil Horizontal Subgrade 
large bend radius (large multiplier, will increase the effective 


1.0 


= 
2 


a 


50 


Bend Radius Multiplier, N 
FIG. and Maximum Pipe Movements 


soil bearing area. Higher coefficient subgrade reaction, will also increase 
the soil resistance bend. Fig. shows the influence both and k,, 
the maximum pipe movement for 16-in. (406-mm) steel pipe with 
500 (152.4 heated temperature differential 140° (60° C). 
Note that the greater the frictional resistance (C-value) between soil and coating, 
the smaller the pipe movement. The greater the soil horizontal subgrade 
reaction coefficient (k,,) the bend, the smaller the pipe movement. 
relationship such this can used assess the cost-effectiveness selection 
trench backfill soil and compaction requirements, because the soil type and 
density influence the and k,, values. addition, for given coating and 
soil type, the pipe movement can estimated increasing bend radius. Thus, 
with proper choice coating and appropriate design and construction, information 
such Fig. may used limit the potential pipeline movements without 
damage the coating. However, for the time being, little known the 
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coating integrity when subjected the thermal induced cyclic movements 
during the design life the pipeline. 

From geotechnical engineering viewpoint, depends soil type, density, 
and depth pipe burial. also varies with soil deformation and strength (2). 
Since heated oil pipelines are characterized low operating frequency (about 
times/yr) cyclic loading acting generally sandy trench backfill 
shallow depths m-3 m), can estimated from nonlinear, drained stress-strain 
relationship sands. The range shown Fig. represents the estimated values 
concrete sands backfill compacted trench about 2.5 deep (8). 

The value coating dependent and also varies with soil type, density, 
and operating temperature. Modified direct shear tests heated samples 


a 
a 


Distance to Bend, d, feet 
FIG. Casing for 1,000 Pipe Two Operating Temperatures 


given coating against sands were conducted (8). The range values (105 
partly from results such tests and partly from 
field average coating wrinkle spacing, (note that emphasized 
that these values are preliminary and indicate only order magnitude 
because the many limitations associated with testing coating using direct 
shear device. PCSI (Pipe-Coating-Soil Interaction) tester being developed 
(9) determine and evaluate pipe coating performances elevated 
temperatures. 

Influence Casing.—Since the pipe within casing not subjected shearing 
resistance like the pipe buried trench, pipe movements near the casing 
entrance well near bends are locations high damage potential. Fig. 
shows this phenomenon two different operating temperatures. Fig. 
represents 16-in. (406-mm) steel pipe buried with 42; 
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known subjected damage when the cyclic movements exceed 0.5 in. 
(13 mm), pipeline designed from Fig. will subjected failure within 
the reach about 100 (30.5 from the bend. Fig. also shows that 


pipe failure can expected the casing entrance the operating temperature 
differential exceeds 180° (82° C). 


Comparison Data 


16-in. (406-mm) steel pipe with 90° bend and but casing 
was instrumented monitor pipe movements (9). The soil cover ranged from 
m-3 the bend area with backfill ranging from concrete sands 
silty sands. Temperature the bend was not exactly known but the temperature 
delivery and receiving stations ranged between 150° and 160° (66° +C- 
71° +C). Fig. shows comparison between calculated pipeline displacements 


d, meters 


|Monitoring Points| 
feet apart 


Distance from Bend, d, feet 


FIG. 11.—Comparison between Pipe Movement Profile and Measured Field 
Data 


and field measured data. Detailed analysis given Ref. Considering many 
uncertainties that involve soil properties, construction procedures, and pipeline 
operations, the comparison shown Fig. considered satisfactory and 
generally substantiate the soil-pipe interaction analysis. 


Based field observations damaged pipe, soil-pipe interaction analysis, 
and field measurements, the following conclusions appear appropriate: 


The expansion and contraction heated oil pipelines produce soil-pipe 
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interface shear and cyclic movements. The magnitude the interface shear 
and movements depends soil properties and pipeline design such length 
straight pipe sections, bend radius, and casing length. 

Excessive interface shear and movements may damage protective coatings 
and expose the pipeline and lead corrosion-induced failures. 

The interface shear and pipe movements can estimated using the soil-pipe 
interaction equations which have two principal parameters. One the soil 
coefficient horizontal subgrade reaction, k,,. The other parameter appears 
unique heated oil pipeline problems. the interface shear parameter, 
and its related relative stiffness parameter Parameter measures the 
coating resistance pipe movement elevated temperatures. The interface 
shear parameter coating dependent and more research needed. 

Increase and values will decrease pipe movements thereby decreasing 
pipeline damage potential. For given pipeline geometry and soil conditions, 
the movements calculated using the soil-pipe interaction analysis showed promise 
one the criteria for selection coatings for heated oil pipelines. 


ACKNOWLEDGMENT 


Portions this work were funded the Southern California Edison Company. 


Their permission publish acknowledged. Allen reviewed the manuscript 
this paper. 


Argent, T., Role Soil Assessment Techniques the Selection Pipeline 
Proceedings the Second International Conference the Internal and 
External Protection Pipes, BRHA Fluid Engineering, 1977. 

Audibert, E., and Nyman, J., Resistance Against Horizontal Motion 
Pipes,’’ Journal Geotechnical Engineering Division, ASCE, Vol. 103, No. GT10, 
Proc. Paper 13303, Oct., 1977, pp. 1119-1143. 

Krizek, J., and McQuade, V., Buried Concrete Journal 
the Geotechnical Engineering Division, ASCE, Vol. 104, No. GT7, Proc. Paper 
13899, July, 1978, pp. 815-837. 

Lange, M., discussion Corrosion Protection Underground Pipelines 
Plastic Tape Pipeline Harris, Proceedings the Second International 
Conference the Internal and External Protection Pipes, 1977, pp. x-13. 

Ligon, B., and Mayer, R., Resistance Buried Pipeline Coatings 
Pipeline and Gas Journal, Feb., 1975, pp. 33-35. 

Strauss, D., and O’Keefe, W., for Nuclear Power, June, 1978, 
pp. 

Wang, L., and O’Rourke, Buried Lifeline Earthquake 
Current State Knowledge Life-Line Earthquake Engineering, ASCE, 
Aug., 1977, pp. 252-266. 


Yen, Interaction Heated Oil Report Southern California 
Edison Company, Dec., 1978, pp. 1-89. 


The following symbols are used this paper: 


effective bearing area pipe bend, L’; 
length casing measured from zero pipe movement, 
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soil-pipe interface shear parameter, 

coefficients for determining pipe movements, 

modulus elasticity the pipe 

coefficient horizontal subgrade reaction, 

pipe length from zero movement pipe bend; 

ratio bend radius pipe radius; 

pipe internal pressure, 

soil reaction force bend pipe expansion direction, 
interface shear resistance per unit length pipe, F/L; 

radius pipe, 

thickness pipe wall, 

pipe movement location 

coefficient thermal expansion, temperature; 
soil-pipe interaction parameter, 

apparent friction angle between soil and coating, degree; 
temperature differential between hot oil and ambient soil, degree; 
unit weight soil, and 

Poisson’s ratio, dimensionless. 
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HEAD PERFORMANCE CONCRETE 
(Reviewed the Pipeline Division) 


The performance limits concrete pipe, specifically prestressed concrete 
cylinder pipe (PCCP), are continually being extended higher and higher internal 
pressure heads. order withstand these heads, pipe cores must heavily 
prestressed during manufacture, and steel spigot rings must carefully analyzed. 
This all requires that the stresses the pipe core and spigot ring accurately 
predicted from manufacture through joining through pressurization. 

Core prestress, which often requires multiple layers wraps prestressing 
wire heavily prestressed pipe, causes the core deflect radially inward. 
This induces initial residual longitudinal bending stresses the steel spigot 
ring. When pipe are joined, the spigot ring compressed radially inward 
the gasket and ring contact the much stiffer bell, also inducing longitudinal 
bending stresses. Under internal pressure, the pipe core and joint expand radially 
but not uniformly, because the difference stiffness the joint area, thus 
causing additional bending stresses. 


The theory presented this paper based the assumption that the joints 
and barrel parts the pipe core can modeled relatively flexible spigot 
attached pipe barrel, restrained from excessive outward radial deflection 
under internal pressure rigid bell attached another pipe barrel. example 


this shown Fig. the bell the right and the spigot the 
left. 


the April 2-6, 1979, ASCE Convention and Exposition, held Boston, 
Mass. 


Research and Development, Hitco, 1600 West 135th St., P.O. Box 1097, Gardena, 
Calif. 90249. 

Note.—Discussion open until June 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication April 23, 1979. 
This paper part the Transportation Engineering Journal ASCE, Proceedings 
the American Society Civil Engineers, ASCE, Vol. 107, No. January, 1981. 
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earlier version the theory was developed Gordon Robertson 
Ameron, and appears proprietary unpublished paper. Tests 108-in. 
(2,740-mm) inside diameter (ID) prestressed concrete cylinder pipe were conduct- 
and was found that the theory needs account for the radial spigot 
ring restraint produced the bell part the joint. The 189-in. (4,800-mm) 
PCCP tests reported the paper verify the expanded theory. computer 
program has been written solve the complex equations multiple locations 
along the total pipe. 

The loads that act produce stress and deformation the spigot result 
from prestressing load, radial gasket load, internal water pressure, external earth 
loads, and radial restraint load acting through the raised flat portion the 
spigot ring that contact with the bell ring. Calculation the latter load 
requires knowledge the pressure induced radial deflections the bell and 
spigot their point contact. 

The method analysis used determine the response the joint components 
the imposed loads based the theory beams elastic foundation 


FIG. 1.—Typical Prestressed Concrete Cylinder Pipe Joint Configuration 


(2). This theory, applied pipe, assumes that longitudinal element 
relatively thin pipe under axially-symmetric loading may regarded 
beam elastic foundation, with the elastic supports along the beam being 
provided circular elements pipe. For simplicity, pipe wall longitudinal 
element unit width analyzed. 

account for discontinuity the pipe wall longitudinal element the bell 
spigot ends, the beam the junction the end and barrel. The 
parts are then analyzed separately for deflections and rotations the cut caused 
known external loading along the beam, and unknown continuity moments 
and shears the cut. equating deformations the cut, the continuity moments 
and shears can determined. Knowing these values and the response external 
loading, final deflections and stresses can calculated the end and the 
barrel. 

The short bell end the pipe, which has embedded steel bell ring and 
nearly thick the barrel, analyzed rigid beam constant thickness. 
The spigot end the pipe analyzed flexible beam, also constant 
thickness, with credit given for the mortar lining. 
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The purpose this paper report tests large diameter PCCP carried 
out verify existing pipe design theory (1) and new unified theory, based 
part elastic analogy, developed for predicting 
pipe core and steel spigot ring stresses from manufacture pressurization. 

Two 189-in. (4,800-mm) lengths PCCP with four layers prestressing 
wire were pressurized 1,386 (422 head 600 psi (4,137 kPa). 
The spigot ring the joint between the two pipes was strain-gaged, were 


each layer prestressing wire both pipe. Strain-gage data showed good 
correlation with predicted values. 


Manuracture 


Overall test pipe design dimensions were: (1) 189 in. (4,800 mm); (2) 
core thickness in. (381 mm); (3) overall length ft, 4-1/2 in. (3.16 
m); (4) spigot length 6-1/4 in. (159 mm); and (5) bell length 5-3/4 in. 
(146 mm). typical PCCP joint shown Fig. 
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FIG. 2.—Strain-Gage Locations Steel Spigot Ring in. 25.4 mm) 


Four layers 1/4-in. (6.35-mm) diam., American Society for Testing and 
Materials, ASTM A648 Class III prestressing wire were applied tensile 
load approx 8,900 (4,036 kg), producing 3,430 psi (23.6 MPa) prestress 
the pipe core. After each layer prestressing wire was applied, strain gages 
were attached the wire and strain readings were taken determine the 
effect applying successive layers wire the previously applied layers. 

The strain-gaged spigot ring had 1/16-in. (1.59-mm) gages along its length 
oriented the longitudinal and circumferential directions, shown Fig. 
They were read before and after each layer wire was applied. 

Gages the prestressing wire and spigot ring were located four points 
around the pipe circumference, shown Fig. 

Typical plots longitudinal and circumferential spigot strains for actual test 
data and values predicted the new theory are shown Fig. Stresses 
can assumed times microstrain. 

Predicted average change prestressing wire strain due additional layer 
wire being applied was 216 microstrain. Test values averaged 201 microstrain, 
which was within the predicted value. 
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few data points Fig. that were not consistent with the others were 
due bad strain-gage readings caused problems maintaining even resistance 
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FIG. 3.—Strain-Gage Locations Around Pipe Circumference in. 25.4 mm; 
304.8 mm) 
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FIG. 4.—Typical Longitudinal Bending and Circumferential Strains Induced Steel 
Spigot Ring from Layer Prestressing Wire in./in. mm/mm) 


the gage attachment the gage meter following each application each 
layer prestressing wire. 
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Pressure Test 


diagram the pressure test setup shown Fig. Pipe with the 
strain-gaged spigot ring, the bottom. 

inside displacement tube, leaving (25-mm) clearance between its 
outside diameter and the pipe’s inside diameter, was used minimize end 
thrust the bulkheads the test setup. Special floating ring bulkhead sealing 
components were used approximately match adjacent pipe joint stiffness 
the top and bottom bulkheads. Details are shown Fig. 


FIG. 6.—Detail Diagram Showing Components Displacement Tube Bulkhead Parts 


Predicted versus actual longitudinal bending strains the spigot after being 
joined the bell, assuming 200-lb/in. (35.7-kg/cm) gasket load are plotted 
Fig. 

Predicted versus actual values tensile strain the prestressing wires versus 
internal pressure plotted for selected gages Fig. These pertain the 
first full pressure cycle. The predicted values are conservative. 

Predicted versus actual longitudinal bending strains the steel spigot ring 


due internal pressure are plotted Fig. for 450 psi (3,100 kPa) 
internal pressure. 


\ 
FIG. 5.—Pressure Test Setup 
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slightly over 450 psi (3,100 kPa), yielding the spigot ring occurs. This 
shown quite dramatically Fig. 10. this figure, the predicted yield strain 
plot not set value because there are residual strains the spigot due 
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FIG. 7.—Predicted Versus Actual Longitudinal Bending Strain Steel Spigot Ring 
Due Joining in./in. mm/mm; in. 25.4 mm) 
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FIG. 8.—Predicted Versus Actual Prestressing Wire Strains Due Pressure in. /in. 


manufacture and joining that are present prior pressure tests. the 600-psi 
(4,137-kPa) pressure level, excessive spigot yielding occurred. This was clearly 
evident after the test pipes were pulled apart and there was seen pronounced 
bulge the spigot ring. 
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FIG. 9.—-Predicted Versus Actual Longitudinal Bending Strains Steel Spigot Ring 
Due Internal Pressure in./in. mm/mm; psi 6.8948 kPa) 
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FIG. 10.—Longitudinal Bending Stresses Spigot Ring Due Internal Pressure 
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Interestingly, succeeding pressurization cycles 600 psi (4,137 kPa), 
shown Fig. 11, indicate that nearly all the spigot ring yielding occurred 
the first cycle. Apparently, this due the fact that once the spigot ring 
yielded out the much stiffer bell ring, the spigot ring had more room 
yield successive pressurization cycles. 

leaks occurred either the pipe wall the joint area the maximum 
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FIG. 11.—Longitudinal Bending Strains Steel Spigot Ring 600 psi (4,137 kPa) 
for Four Pressure Cycles in. 25.4 mm; in./in. mm/mm) 


test pressure 600 psi (4,137 kPa). Subsequent mortar grouting the test 
pipe joint with mortar placed the space between spigot ring and bell ring 
showed further reduction spigot ring stresses during pressurization. 


good correlation between predicted and measured strains the prestressing 
wire and the steel spigot ring due prestressing, joining, and pressurization 
appears verify the new unified theory developed. Thus, the new theory may 
safely used the design alternate spigot ring configurations. 

noted that although current design practice limits steel spigot ring 
longitudinal bending stresses the yield strength the steel, such limitation 
not necessary. After spigot ring yields and contacts the bell ring, more 
yielding can take place, and the joint should even tighter because yielding. 
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PERMEABILITY LIME STABILIZED SOILS 


(Reviewed the Highway Division) 


INTRODUCTION 


Many soil types can stabilized economically with lime. Soil-lime mixtures 
may used sub-hase base course for roads, airfields, parking lots, 
highway shoulders, and similar applications (1,2). Soil stabilization with lime 
has good prospect Iraq limestone deposits are available number 
areas throughout the country and the weather suitable for this type 
stabilization (3,4). 

The improvements the engineering properties soil lime added 
are attributed two basic reactions (2): 


Colloidal-type reactions which occur almost immediately and result 
immediate changes soil plasticity, workability, and swell properties (5). 
involves one the following: (1) Ion exchange calcium for the ion carried 
the soil; (2) depression the double layer the soil colloids result 
the increased cation concentration the pore water; and (3) expansion 
the double layer the soil colloids from the high lime. 

Pozzolanic reactions between the lime with the available reactive alumina 
silica the soil, both. 


When lime added fine grained soil, generally increases the plasticity 
index low plasticity soils and decreases the plasticity cohesive soils. They 
become friable and easy work (6,7). Lime reduces the maximum dry density 
and increases the optimum moisture content clayey soils given compactive 
effort (4,7,8). The swelling characteristic the soil reduced and the shrinkage 
limit increased when lime added (6). The strength properties Iraqi soils 
are improved when stabilized with lime (9,10,11). 

Prof., Dept. Civ. Engrg., Coll. Engrg., Univ. Baghdad, Baghdad, Iraq. 

Asst. Instr., Inst. Tech., Baghdad, Iraq. 

Note.—Discussion open until June 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication April 1980. 
This paper part the Transportation Engineering Journal ASCE, Proceedings 
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measure the ease with which fluid passes 
through material. one the fundamental engineering properties 
soil. enters into almost all seepage, settlement, and stability calculations. 

The coefficient permeability used indicate the permeability soils 
(12,13). The coefficient permeability for given soil depends particle 
diameter, properties the pore fluid, the void ratio, the shape and arrangements 
the pores and the soil particles, and the amount undisolved gas 
the pore water (14,15,16,17). 

The permeability compacted soils varies with the molding water content. 
dry side optimum moisture content the permeability high due the 
flocculant soil structure and the low density. marked decrease permeability 
accompanies the increase molding water content the dry side the optimum. 
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FIG. 1.—Grain Size Distribution 


The permeability reaches minimum slightly above the optimum moisture 
content, after which slight increase permeability occurs. The large decrease 
permeability with increasing molding water content the dry side the 
optimum caused the orientation the clay-size particles. Wet optimum, 
the permeability increases slightly because the effects decreasing dry density 
more than offset the effects the particle orientation (2,14,16,18). 

When lime added clayey soil, the soil becomes flocculated and the 
permeability increases. granular soils, the lime reduces the void ratio and 
the permeability. 

The purpose this paper investigate the validity this hypothesis. 
The effect lime the permeability two compacted Iraqi soils was 
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investigated. The effects such factors lime content, molding moisture 
content, and curing age the coefficient permeability were investigated. 


Soils.—Two soils were used, i.e., poorly graded fine river sand from Baghdad 
which typical Iraqi river sands, and sandy silty clay from Tikreet. Fig. 
shows the grain size distribution the two soils. Some the properties 


TABLE 1.—Soil Properties 


Sample Number 


Variable 


(1) 


Location 
Unit weight solids 
Atterberg limits 
liquid limit, percent 
plastic iimit, percent 
plasticity index, percent 
Soil description 


N.P. 

poorly graded, 
gray, fine, river 

sand 


reddish-brown 
sandy silty clay 


TABLE 2.—Chemical Analysis Lime 


weight, percent 


CaO 
MgO 

SiO, 

Na,O 

Loss ignition 


the soils are given Table 
Lime.—Hydrated lime manufactured from limestone deposits from Heet area 
was used throughout this investigation. Table gives the chemical analysis 
the lime used. 
Water.—Ordinary drinking water was used mold specimens various 
moisture contents. Distilled water was used for the permeability tests. 
Compaction.—The standard Proctor compaction hammer and mold was used 
prepare in. (10.1 cm) diam 4.6 in. (11.7 cm) high specimens the 
river sand. The soil was compacted three layers with blows/layer. The 
height fall was in. (30 cm) using 5.5 (2.5 kg) drop hammer. 
Smaller specimens with 2.8 in. (7.1 cm) diam and 1.1 in. (2.8 cm) high using 


(3) 
Baghdad Tikreet 
2.66 2.69 

Oxide 
(1) (2) 
64.96 
7.58 
0.71 
0.64 
0.53 
0.32 
0.16 
24.70 
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model hammer 1.8 (0.82 kg) with drop height in. (20 cm), were prepared 
the sandy silty clay. The compactive effort was equivalent that standard 
Proctor compaction test. was observed that the maximum dry density using 
the model hammer was less than that the standard Proctor test. The 
optimum moisture content the raw soil was about the same shown 
Fig. 

Specimen Preparation.—The required amount lime (percent the dry weight 
the soil) was mixed with the oven dried pulverized soil. Water for the compaction 
was than added. After compaction the specimens and the molds were wrapped 
cellophane paper and sealed with paraffin wax and stored curing room 
100% relative humidity. After the curring and after the wax had been removed, 
the permeability was determined using the compaction cylinder part the 
permeameter reduce the disturbance along the edges the specimens 
they are removed from the compaction cylinder. 
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FIG. 2.—Effect Lime Content Compaction Properties Soil from Tikreet 


Permeability compaction cylinder was used permeameter 
mentioned before. falling head permeability test was used throughout this 
investigation. The clayey soil specimens were subjected pressure 
kg/cm’ for before the start the tests saturate the specimens 
(19). The tests were carried out according British Standards No. 1377-1975 
and room temperature 25°C. The values reported were the average 
three tests. 


River Sand.—The maximum dry density the compacted river sand increased 
with increasing lime content shown Fig. 3(a). This was expected since 
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the river sand was uniformly poorly graded. 


The permeability tests the samples which were investigated just after the 
addition lime were carried out after the molding. This time was needed 
saturate the samples. 


Immediate Tests.—The following results were obtained from the tests carried 
out after molding: 


The permeability was relatively high the dry side the optimum water 
content. decreased with increasing molding water content reaching minimum 
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FIG. 3.—Effect Molding Water Content and Lime Content Permeability River 
Sand 


value slightly wet optimum, after which the permeability increased shown 
Fig. 

The permeability decreased with the increasing lime content for all investi- 
gated water contents shown Fig. 3(b) and The lime acted filler 
plugging the voids and reduced the flow channels the mixture. 

The specimens molded the dry side the optimum had higher 
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FIG. 4.—Effect Lime Content Permeability River Sand Tested 25°C 
after Molding 
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permeability than the specimens prepared wet optimum indicated Fig. 
There certain molding water content that produces minimum permeability 
for all investigated lime contents. This water content was 19%, which 
higher than the optimum water content for this soil-lime mixture. 


Cured Specimens.—For the soil-lime specimens which were cured 100% 
humidity and 25°C, the following results were obtained: 


The permeability was reduced with curing age all molding water contents 
shown Fig. One hypothesis explain this that when the lime hydrates 
with age reacts with the silica the soil, and the gel will grow with curing 
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FIG. 6.—Effect Molding Water Content and Lime Content Permeability Tikreet 
Soil (One Day after Molding) Tested 25°C 


age filling the voids. The decrease permeability could also due the 
precipitation calcium carbonate the voids the ionized calcium reacts 
with the dissolved carbon dioxide the water. 

There exists molding water content which the permeability minimum 
regardless curing age. This molding water content about higher than 
the optimum moisture content shown Fig. 

Sandy Silty Clay.— When lime was added the clayey soil, the maximum 
dry density decreased and the optimum moisture content increased with increas- 
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ing lime content shown Fig. This agrees with other results reported 
the literature (2,4,9). The permeability tests the compacted clayey soil- 
lime mixture gave the following results: 


The permeability the compacted soil without lime decreased with 
increasing molding water content. reached minimum water content 
slightly the wet side the optimum moisture content. The permeability 
increased when the water content was increased further shown Fig. 6(a). 
This agrees well with the results reported from other tests the literature. 
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FIG. 7.—Relations between Dry Density Soil from Tikreet Stabilized with Lime 
and Permeability (Tested One Day after Molding 25°C) 


The density was low dry optimum and the soil flocculated, thus the 
permeability high. the molding moisture content increased the density 
increased and the soil particles become more oriented resulting less permeable 
soil. When the moisture content increased further (wet optimum), although 
the soil particles are well oriented (dispersed structure), the permeability increases 
shown Fig. 6(a). This due the low density wet optimum. 

When lime was added, the soil became flocculated. When gel formed around 
the soil particles, the soil particles grew size and the voids the mixture 
increased compared with raw soil. Therefore the permeability the soil-lime 
mixture was higher than that the original clayey soil without lime shown 
Fig. 6(a). 
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and lime there was molding water content which the 
permeability reached maximum. This molding water content was identical 
the optimum moisture content for that soil-lime mixtures shown Figs. 

Fig. 6(b) indicates that lime content did not affect much the permeability 
soil-lime mixtures dry optimum. This probably because the soil already 
flocculated structure. The largest effect the lime permeability occurred 
the optimum moisture content. The coefficient permeability specimens 
one day after molding increased with the addition lime where 
the permeability reached maximum, after which the permeability decreased 
lime. increased again lime. The effect lime the compressive 
strength showed similar pattern with maximum lime (10,11). Specimens 
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which were compacted wet optimum moisture content behaved similar 
manner that compacted optimum. 

For compacted soil without lime the coefficient permeability decreased 
with increasing dry density should expected. When lime was added, the 
permeability increased with increasing dry density. Fig. indicates linear 
trend 0%, 5%, and lime. the same dry density the permeability increased 
with increasing lime content, probably due the increased flocculation higher 
lime contents. 

The permeability soil-lime mixture decreased with curing age for all 
moisture contents tried shown Fig. 
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At3% lime the permeability reached maximum slightly wet the optimum 
molding water content for all curing ages shown Fig. 


The following conclusions may drawn, limited the soils and test conditions 
investigated: 


The permeability characteristics soil-lime mixture depends the type 
soil. For granular soils the permeability decreases with the addition lime 
while for cohesive soils the permeability increased. 

The molding water content soil-lime mixture influences the coefficient 
permeability. How the coefficient permeability varies with the molding 
water content depends the type soil being granular fine grained. 

For granular soil-lime mixtures there molding water content which 
the permeability minimum regardless the lime content. This molding 
water content was found slightly the wet side the optimum. 

For cohesive soil-lime mixture there molding water content that produces 
maximum permeability. This corresponds about the optimum water content. 

For granular soils the permeability decreases with increasing lime content 
all molding water contents investigated. 

For both soil types the permeability decreased with curing age whether 
the mixture was prepared dry optimum, optimum, wet optimum. 

For cohesive soil-lime mixture the coefficient permeability increased 
with increasing dry density. 

same dry density, the permeability increased with increasing lime content. 
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Tax ALLOCATION FORMULAS 
STEPWISE REGRESSION 


Lansford Bell," ASCE and Russell Shelton” 
(Reviewed the Highway Division) 


INTRODUCTION 


significant source operating capital for city and county road programs 
statewide gasoline motor fuel tax revenues. portion the statewide 
revenue typically distributed county, city, other local municipality 
according its population, road mileage, surface area, traffic volume, some 
combination these other variables. 

1979, legislation was introduced states increase statewide gasoline 
tax (7). Some the legislation introduced 1979 contained provisions for 
altering the distribution formulas that provide revenue local road programs. 
bill containing provisions alter Alabama’s distribution formula was considered 
and then rejected committee action during the 1979 legislative session. Similar 
proposals number other states have been reported the writers. 

Revenue distribution formulas have evolved over the years the result 
political compromise with little attempt rational analysis. Many distribu- 
tion formulas are outdated that they were established time when there 
was pressing need construct new roads rural areas facilitate transporting 
agricultural products into urban market centers. Today, however, great 
percentage local road expenditures toward roadway maintenance, bridge 
replacement and repair, and safety improvements. 

previous research project conducted Auburn University (3) compiled 
data the current condition Alabama’s secondary county road system. 
The research results cited need for $0.04/gal increase the statewide 
gasoline tax, assuming the counties were continue receiving 55% all collected 
revenue. became apparent during the course the research, however, that 
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revenue generated the tax increase would not coincide with the critical needs 
the counties the additional revenue was distributed according formulas 
being considered the state legislature that time. 

second research project was therefore initiated develop formula that 
best allocated revenue generated increase statewide gasoline tax 
Alabama’s county road departments. Six variables thought predictors 
county revenue need were initially selected for possible inclusion the 
formula. Stepwise regression techniques were then applied the data select 
the optimum combination variables and their respective percentages the 
formula. Presented herein examination the procedure used determine 
proposed allocation formulas and assess their adequacy. 


Specific allocation formulas are used states allocate motor fuel tax 
revenue county road programs. Specific formulas are used states 


TABLE 1.—Variables Used Distribution Statewide Motor Fuel Receipts 


Number states 
Number using variable 
Variable used states using part 
determine variable part distribution formula 
statewide distribution formula counties and municipalities 


(1) (2) 


Road miles 

Population 

Vehicle registration 
license fees 

Land area 

Revenue 

Need 

Vehicle miles 

Property values 

Farm products sold 

Linear feet bridges 


allocate motor fuel tax revenue cities other municipalities. The most common 
variables used motor fuel allocation formulas are shown Table (9). 
should noted that most states combine two more variables their formula 
and many states allocate some portion revenue equally all counties 
municipalities. 

workable tax allocation formula should ideally possess the following 
characteristics: 


Variable data should readily available inexpensive tabulate. 

Variable data should not subject manipulation the local jurisdictions. 

The formula shouid not unnecessarily complex. Linear combinations 
variables should employed unless can shown that more complex 


TE1 TAX ALLOCATION FORMULAS 


expressions result higher correlation with revenue-need. 
The formula should compatible with the purpose for which the tax 
tax increase imposed. 


the writers’ knowledge, the only comprehensive investigation motor 
fuel allocation formulas was performed the state Washington (5,10,11,12). 
This investigation was conducted the University Washington and the State 
College Washington using funds provided the Washington legislature. 
The research resulted formula that distributed revenue the state’s 
counties follows: 10% equally; 30% the basis vehicle registration; 
30% the basis trunk mileage multiplied the prorated estimated annual 
cost per trunk mile; 30% the basis trunk mileage multiplied prorated 
money need factor. The estimated cost per trunk mile and the money need 
factor are updated periodically. 

The state Washington has apparently achieved the desired level coordina- 
tion between university researchers and political advisory committee. Recent 
correspondence from Washington highway officials indicates, however, that 
changes achieve greater degree formula simplicity are being contemplated. 


Method.—Multiple linear regression analysis can used predict, estimate 
the value one dependent variable from known values other independent 


variables. The method involves the assumption linear relationship the 
form 


which linear function predictor variables X,, X,, ..., X,. The 
coefficients and the error term, are never known exactly. The researcher 
therefore postulates prediction regression equation the form: 


which the predicted value and the regression coefficients and 
estimates and The overall goal regression analysis obtain 
values for the regression coefficients and then perform statistical tests evaluate 
the adequacy the model. 

Stagewise regression procedure which large set predictor variables 
can narrowed down eliminating variables that not make significant 
contributions the regression equation. Forward elimination, backward elimina- 
tion, and stepwise regression are techniques for successively examining combina- 
tions predictor variables and then revising the variable set according 
predetermined statistical tests. 

Fuel Tax Allocation Model.—Consider the case where highway revenue need 
can accurately determined estimated for counties political jurisdictions 
within state. Assume data defining population, road mileage and other possible 
need predictors are also available for each county. Regression techniques could 
then used determine which variables what relative proportion should 


included linear equation best distribute fuel tax revenue throughout 
the state. 
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obvious question asked not distribute revenue according 
need data reflecting need can The answer simply that 
formulas involving predictor variables such those shown Table are simple, 
acceptable the public and not subject political manipulation debate. 
Revenue need is, however, somewhat subjective and may function 
local priorities. The need variable employed the regression analysis should 
ideally calculated using weighting factors, priorities guidelines approved 
state legislative committee, county associations other designated organiza- 
tions. 

Computer Software.—Multiple regression analysis computer programs are 
available most large computing facilities. The BMDP Biomedical Computer 
Programs (4), the SPSS Statistical Package for the Social Sciences (6), and 
the SAS Statistical Analysis System (2) are each capable performing regression 
analysis and stagewise variable selection. The use these program packages 
and their relative merits are also examined Ref. 13. 


Model Variables.—Alabama’s current fuel tax distribution formula, enacted 
1967, allocates 55% all statewide $0.07/gal gasoline tax revenue the 


TABLE 2.—Variables Used Regression Analysis 


Variable Description 
(1) (2) 

POPL County population, 1978 estimate 

MILE Miles county roads, 1978 

AREA Square miles county land area 

VEHR Total number vehicles registered county, 1978 

TRUCKR Total number trucks registered county, 1978 

MTLS Thousands tons timber, coal and other bulk materials 
mined county, 1978 

NEED $13,000 times miles road need immediate resurfac- 

(dependent variable) ing plus cost associated with all bridges need imme- 

diate replacement 


counties such that 45% the counties’ share distributed equally and 
55% distributed the basis population. Previous research indicated that 
legislation proposing increase statewide gasoline tax should include 
improved allocation formula. 

first step deriving distribution formula for the additional revenue, 
six predictor variables, defined Table were identified. The variables POPL, 
MILE, AREA and VEHR are the variables most commonly used other states. 
The variables TRUCKR and MTLS were included because local haulers 
timber, coal and other bulk materials have apparently accelerated road deteriora- 
tion many counties. Data for each variable were then compiled from various 
state agencies. 

The dependent variable, NEED, was determined from data compiled 
1979 questionnaire survey county engineers conducted part the research 
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described Ref. was felt that any increase tax revenue would 
spent primarily resurfacing and bridge replacement. The value the need 
variable was therefore computed $13,000 times the number miles road 
need immediate resurfacing plus the cost associated with all bridges 
need immediate replacement. 

Defining need this manner results $307,000,000 total revenue requirement 
for all counties. This $307,000,000 estimated revenue requirement would 
generated about the current $0.07/gal gasoline tax was raised 

Regression Analysis.—The regression analysis model can now stated 
equation estimate need function six predictor variables: 


which NEED the estimated revenue requirement for given county; and 
the coefficients b,, b,, ... are determined regression analysis. 

obtain equation terms only those variables that contribute 
significantly the prediction revenue need, variable data were entered into 
the BMDPIR stepwise regression computer program. Variables are entered 
removed performing partial tests which determine should 
rejected some stage the variable selection procedure. When cannot 
rejected giving the same information some other predictor variable 
the formula. The level significance which variable entered removed 
determined the program user. significance level 0.05 was used 
herein. For additional information relating statistical tests, variable selection 
and program controls the reader referred Ref. and 13. 

initial program execution produced the regression equation: 


which NEED the predicted allocation given county units dollars. 
Removing VEHR from the data set attempt obtain formula 
terms the five remaining variables results the regression equation: 


NEED $750,180 125(TRUCKR) 2,960 (MILE) 


And finally, removing TRUCKR from the data set attempt formulate 
equation terms POPL, MILE, AREA, and MTLS results in: 


NEED $415,420 14.10(POPL 3,840 (MILE) 


interesting note that MTLS and AREA were not selected for entry 
into Eqs. 4-6. There therefore 95% confidence that these variables 
not contribute significantly the estimation prediction need. 

Formula Adequacy.—Two standard statistical parameters can calculated 
assess formula adequacy compare proposed allocation formula with 
one current use. The standard error the estimate measures the amount 
that actual need varies from need predicted the regression equation. Standard 
errors associated with Eqs. 4-6 were calculated $2,430,000, $2,440,000, and 
$2,450,000, respectively. These values can interpreted mean that 
$307,000,000 was distributed counties using Eq. then the allocation 


7. 
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would vary from actual need $2,430,000 for average county. 

The standard error the estimate can used compare regression formula 
formula current use other formulas under consideration state 
legislature. Regression Eq. compared two such formulas Table 
can seen that the standard error the regression equation significantly 
less than the standard error Alabama’s current allocation formula. 

second statistic that assesses the overall adequacy regression equation 
the multiple coefficient determination, Values associated 
with Eqs. 4-6 were within range 0.43-0.45 indicating that only 43%-45% 
the variation need can explained the regression equation. This could 
interpreted justification for distributing portion fuel tax revenues 
according documented need when can accomplished without political 
bias data manipulation. 


TABLE 3.—Standard Error Alternate Allocation Formulas 


Variables 
Allocation (relative weights), Standard error, 
formula percentage millions dollars 
(1) (2) (3) 
Regression equation (equal distribution) 
(by population) 
(by road mileage) 
Traditional 1/3 formula 33(by population) 
(considered previous (by road mileage) 
Alabama legislation— (by land area) 
not enacted) 
Alabama’s current formula (equal) 
(by population) 


Eqs. 4-6 employ different variable combinations but produce standard test 


Statistics that are approximately equal because the variables TRUCKR, VEHR, 
and POPL are highly correlated. 


Statewide fuel tax revenues are typically distributed local jurisdictions 
the basis population, road miles, registered vehicles, and other variables 
for which data are readily available. Many states are considering legislation 
increase statewide fuel taxes but little attention has been given 
the adequacy current proposed allocation formulas. 

Multiple linear regression analysis well suited the derivation simple 
fuel tax allocation formulas that best distribute statewide revenue local political 
jurisdictions. Although arbitrarily defined the analysis performed herein, the 
dependent variable NEED should ideally structured with input from legislative 
advisory committees county associations. 

Six variables that were thought related current county highway needs 
Alabama were identified and regression analysis performed. Three gasoline 
tax allocation formulas were obtained that would appropriate for distributing 
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revenue generated from proposed increase Alabama’s statewide gasoline 
tax. 


The degree improvement associated with allocation formula derived 
regression techniques best measured the standard error the estimate. 
For total allocation $307,000,000 Alabama counties and the assumed 
definition need, typical regression equation improves the error allocation 
average county ($3,380,000-$2,450,000) $930,000. 
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dependent variable linear model; 
dependent variable regression equation; 
constant linear model; 

coefficient linear model; and 

error term linear model. 


Subscripts 
number predictor variables used regression equation. 
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INTERACTION AND PIPELINE 


(Reviewed the Pipeline Division) 


The sizing straight segments buried pipelines based conventional 
analysis methods. The line first sized react internal pressure using 
the well-known relationship for hoop stresses: 


which o,, hoop stress, pounds per square inch kips per pascal; 
internal pressure, pounds per square inch kips per pascal; average 
pipe radius, inches millimeters; and pipe wall thickness, inches 
millimeters. the line subjected temperature differential, check 
must also made insure that longitudinal stresses and principal shear stresses 
not exceed design allowables. 

Segments buried pipeline some distance from bends the line are fully 
restrained due soil friction. The longitudinal stress the area is: 


which longitudinal stress, pounds per square inch kips per pascal; 
Young’s modulus, pounds per square inch kips per pascal; 
coefficient thermal expansion, inches per inch per degrees Fahrenheit 
millimeters per millimeter per degrees Celsius; temperature differential, 
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degrees Fahrenheit degrees Celsius; and Poisson’s ratio. 
For designs with high temperature differentials, the critical stress usually 
the principal shear stress which obtained from the relationship: 


max 2 

locations where the line changes direction, thrust force developed 
which may result undesirable movement adequate restraint not provided. 
Undesirable movements the pipe are translated into pipe distortions and 
associated critical stresses. The logical approach the design buried pipeline, 
therefore, provide adequate restraint. 

The types restraints used react thrust forces include anchor blocks, 
thrust blocks, soil uplift resistance, bearing capacity, and passive resistance. 
This paper concerned with the study soil the restraining force pipe 
bends. The magnitude movement which can tolerated bend outside 
the scope this paper. 

related papers, Spangler (7) presents treatise the design buried 
flexible conduits with lateral movement caused deformation the flexible 
conduit being restrained backfill material. concludes that the Rankine 
theory not adequate describe the response the behind flexible 
pipes since predicts passive pressure and not the movement required mobilize 
these pressures. His empirical formulation (lowa formula) function the 
pipe geometry, material properties, and the horizontal subgrade reaction the 
backfill. 

The modulus subgrade reaction complex soil property which difficult 
determine with the degree accuracy required obtain reliable results 
with Spangler’s formulation. Audibert and came similar conclusion 
their scale model tests 6-in. pipes. 

Luscher, al. (5) presented paper the pipe soil interaction the 
Trans-Alaska pipeline. Their study pipe-soil interaction based modified 
Rankine equations. 

The study presented herein evaluates the pipe-soil interaction for various 
types backfill and predicts backfill failure pipe bends for rigid steel pines 
in. (1,067 mm) diam. Finite element modeling techniques were used 
the basis for the analysis. The unique advantage this method that basic 
Stress-strain parameters the soil backfill are employed and derived soil 
property (such subgrade reaction) which generally the source error 
and controversy, needed the analysis. Practical design suggestions are 
offered considering various pipe sizes and forces, depth cover, material 
availability, construction economy, and expediency. 

The finite-element method utilized this study yields displacement confi- 
guration for given set design conditions. Knowing the displacement confi- 


guration, movements and stresses can computed establish critical movements 
bend. 


STATEMENT 


free-body diagram fully restrained bend segment pipeline shown 
Fig. The resultant force (F) equal to: 
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which resultant end force; internal pressure, pounds per square 
inch kips per pascal; internal area pipe, square inches square 
millimeters; and steel cross-sectional area, square inches square 
millimeters. 


The pipe soil interaction force, P,, required fully restrain bend segment 
given by: 


which the pipe-soil interaction force, pound feet per foot newtons; 
resultant end force, pound feet per foot newtons; and pipe 
bend radius, feet meters. 


The forces resisting movement for different type bends are given follows: 


Type bend Restraining force 


Horizontal Soil passive resistance 
Sagbend Soil bearing 
Overbend Weight full pipe, 
weight soil overburden, 
and soil friction. 


The development Eq. based the assumption movement 
the bend segment. With the exception the weight full pipe and soil 


FIG. Body Diagram Fully Restrained Pipe Bend 


overburden, movement must occur mobilize the soil restraining forces. This 
movement results the development shears and moments bend segment 
shown Fig. should noted that pipe bend displaced, the 
resultant end force, will decrease magnitude caused relaxation the 
resistance thermal expansion. This reduction force normally neglected 
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since the moments and shears associated with movement large enough 
significantly decrease the resultant end force would yield critical pipe stresses. 

From Fig. apparent that the resistance the thrust force bend 
reacted combination the structural strength the pipe well 
the backfill soil resistance. Since the mobilization soil resistance depends 
pipe movement, insufficient lateral movement very rigid pipe may not 
take advantage the available backfill resistance. the other hand, excessive 
lateral movements may cause failure the backfill and the eventual separation 
the pipe from the surrounding soil. Continued use the pipeline under 
these circumstances may overstress the pipe bend structurally. 


FIG. Body Diagram Pipe Bend with Movement 


This study makes reference horizontal bends but the conclusions apply 
equally overbends and sagbends. 


finite-element approach was used analyze the pipe-soil interaction. The 
computer program described Radhakrishnan and Jones (6) was adapted 
the problem the following manner. 


was recognized that the actual radial force shall vary along the length 
the bend three dimensional manner. However, two-dimensional analysis 
long section the bend subjected maximum radial force per foot 
was considered satisfactory. The weight steel and the weight oil the 
long section the pipe bend were combined and divided the total 
volume the long pipe section arrive density material which 
was mathematical equivalent the previously mentioned two components. 
Other parameters were assigned this mathematically equivalent ficticious 
material consistent with the nonlinear model described the following paragraphs. 
The cross section the pipe was idealized number straight line segments 
since the program does not incorporate curved lines. The pipe-soil interaction 
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force (P,) calculated was distributed equally among the nodes the finite-element 
mesh within the pipe sections. typical finite-element mesh shown Fig. 
the arrows identify the nodes carrying the pipe soil interaction force, P,. 

The soil enveloping the pipe bend was modeled continuum one 
material assigning parameters representative their type determined 
laboratory testing. The nonlinearity the stress-strain relationship was modeled 


hyperbolic transformation suggested Kondner and Zelasko (4) and 
defined the following relation: 


which strain; and major and minor stresses respectively; 
the failure ratio stress difference failure the ultimate asymptotic 
value the hyperbolic stress strain curve; and and the intercept and 
slope respectively the straight line represented the first part Eq. 


RANKINE FAILURE LIMIT 
OF ACTIVE WEDGE 


FIG. 3.—Finite Element Mesh 42-in. (1,067 mm) Pipe Bend 


The confining pressure dependency the stress-strain modulus was modeled 


line with Janbu (3), and Duncan and Chang’s (2) research stated the 
following relation: 


which initial, unloading, reloading modulus; primary loading, 
unloading, reloading modulus number; atmospheric pressure; 
minor principal stress; and the modulus exponent. All the parameters used 
the analysis are summarized Table 

The state stress the finite-element mesh was studied utilizing 
diagram. The parameters and are functions effective principal stresses 


The boundaries the soil continuum were established suf! ficient distance 
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away from the pipe such that the influence the radial force the soil 
continuum was not affected them. relatively small grid (elements) was 
used the vicinity the pipe closely determine the effect movement 


the backfill nearest the source the radial force. 


TABLE 1.—Soil Continuum (Backfill) Parameters 


Variable Dense sand Medium dense sand Loose sand 
(1) (2) (3) (4) 


Unit weight 120 

Poisson’s ratio before failure 0.30 

Poisson’s ratio near failure 0.49 

Coefficient lateral earth 

pressure rest 0.38 

Friction angle 38° 

Cohesion 

2,000 

K unloading 2, 120 
0.55 
0.87 


NOTES 
Friction Angle (Tan 2) 
h / i) 
where, 
3, = major principal stress (effective 
6, = minor principal stress (effective 


FAILURE 


in pounds per square foot (ki lopascals) 


1500 
(23.95) (71.85) 


5, im pounds per square foot (kilopascals) 


FIG. 4.—Soil Backfill Failure Envelopes 


The computer program generates and solves several simultaneous 
for each grid mesh element under the selected nonlinear constitutive stress- 
strain relationship order calculate stresses and displacements the soil 


2000, psf 
(95.80), kPa 


equations 


0.30 0.30 
0.49 0.49 
0.43 0.50 
35° 30° 
720 295 
900 1,090 
0.50 0.65 
0.80 0.90 
g/cm’. 
g/cm’. 
psf, 
kPa (71.85 
(47.90) 
a=26.6° 
STABLE ZOME 
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backfill continuum. The stresses and displacements obtained for various pipe 
sizes, radial forces, soil cover, and backfill types were studied derive the 
conclusions presented the following sections. 


Errect 


The effect backfill the outward radial movement illustrated 
considering three states sand denseness, namely dense sand, medium dense 
sand, and loose sand. Their associated properties are given Table and 
their failure envelope determined laboratory tests are depicted Fig. 


1 Ft. (0.305 m) 
Displacement Scale 


INITIAL POSITION 
OF PIPE 


—~-————. DISPLACED POSITION IN DENSE 


-01 -02 Ft. 
(.0305 m) (.0610 m) 
——--——— DISPLACED POSITION IN LOOSE SAND 


FIG. 5.—Displacement Comparison 42-in. (1,067 mm) Pipe Bend 


42-in. (1,067 mm) diam steel pipe bend 120 (36.6 radius and 
(1.83 soil cover carrying crude oil design pressure 600 psi 
(4,134 kPa) was considered develop outward radial force 8,000 
(116,721 N/m) the bend. Two types soil backfill were studied, namely 
dense sand and loose sand, the resulting outward movement the pipe 
bend depicted Fig. where the displacement scale magnified five times 
relation the geometric scale. The maximum radial displacements were 
calculated 0.09 in. (2.29 mm) and 0.93 in. (23.62 mm) for the dense sand 
and loose sand backfill continuums, respectively. 

The displacements are directly related the stresses any soil continuum. 
general, the stresses the soil backfill elements the inside the bend 
radius will reduce (active case) from the rest condition (no radial force 
development), while the stresses the backfill elements the outside the 
bend will increase (passive case). 

The state stress the finite-element mesh elements exhibiting passive 
case shown Fig. and Fig. for the dense and loose sand continuums, 
respectively, for single load cycle. contrast with the dense sand continuum, 
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the state stress number elements the region significant passive 
pressure plots above the failure envelope for loose sand continuum for the 
same radial force, thus indicating unstable conditions. The maximum radial 
displacement for loose sand was calculated more than ten times that 


920 
foot 


FIG. 6.—Stresses Dense Sand Continuum after Radial Load Application 42-in. 
(1,067 mm) Pipe Bend 


STABLE 
ZONE 


Bor 
16) (32 32) 57. 


FIG. 7.—Stresses Loose Sand Continuum After Radial Load Application 42-in. 
(1,067 mm) Pipe Bend 
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the dense sand. The conclusions about the backfill instability depend 
the backfill type alone and are independent the pipe-bend diameter. 

was observed that the region significant stress increase was 1-2 diam 
front and below the pipe for both dense and loose sand continuums. With 
reference Fig. where Rankine passive wedge extends diam below the 
pipe bottom and intersects the ground surface approximately diam away 
from the pipe, evident that the region significant passive pressure 
development overestimated when compared with finite-element procedures. 
the passive wedge slope were drawn such that the wedge intersects 
the bottom the pipe, i.e., intersecting the ground approximately diam 
away, portion the passive pressure region will underestimated. 
The conventional earth pressure theories not lend themselves appropriately 
defining the region significant stress increase. finite-element formulation 
the other hand appropriately defines the critical area and will eventually 
result better and more economical design. 

The separation pipe from the backfill due high displacements and backfill 
failure condition represents totally new set boundary conditions for the 
stress analysis the pipe. Under these conditions, successive cycles may result 
fatigue failure especially for highly stressed components. Recognition this 
possibility the time design will alert the designer use alternatives such 
better backfill material, anchor blocks, etc. 


Errect Cover 


Overburden pressure function soil cover [Fig. p-g diagram 
(Fig. 4), the initial state stress for deep soil cover will plot the right 
the corresponding plot for shallow soil cover. 

When the radial force applied, given stress path any given point 
the backfill has shorter distance travel before the failure envelope reached 
for shallow soil cover. The mode backfill behavior for shallow and deep 
soil covers schematically represented Figs. 8(b) and 

Analyses were performed for two cases 16-in. (406 mm) pipe with 60-ft 
(18.3-m) bend radius under different soil cover and density conditions. The 
results were: Case I—3 (0.92 soil cover, dense sand backfill, and 
radial force 3,700 (53,984 N/m); and Case II—10 (3.05 soil 
cover, medium dense sand backfill, and radial force 4,000 (58,361 
N/m). Fig. shows comparison displacements for the aforementioned 
two cases during the first cycle loading; the displacement scale magnified 
times illustrate the displacement mode. The zone significant passive 
pressure development was found extend approximately diam-3 diam 
front and under the pipe for both cases. extended the surface for 
Case and diam above the pipe for Case II. The results the analysis 
showed that even though radial force was approximately equal, the state 
stress relatively large number critical elements the dense sand continuum 
(which stronger the two continuums) plotted above the failure envelope 
indicated widespread failure region compared with the medium dense sand 
where only few isolated elements local failure were observed. the writers’ 
opinions, the comparatively larger displacements and the unstable conditions 
the dense sand backfill are the result the shallow soil cover. corollary, 


replacement 


lift out the ground. 


Determination the optimum thickness soil cover and the suitable backfill 


type necessary make the design cost effective. 


inuum Is no 


(b) Shallow Soil Cover; and (c) 
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FIG. 8.—Schematic 
Deep Soil Cover 
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From practical application standpoint, the pipe-soil interaction force (P,) 
should computed using Eq. and compared with the appropriate soil resistance 
available for given subsurface conditions. the available passive resistance 
not adequate, the designer should consider more sophisticated finite-element 
analysis compute the pipe bend movement and the resulting stresses (Fig. 
for proper design. Alternatively, appropriate combination backfill 
type and soil cover can achieved the finite-element analysis considered 
this paper allow for very small bend movement and design conditions 
approaching fully restrained condition (Fig. 1). 


Loading 
Direction 


LEGEND 


DISPLACED POSITION MEDIUM 

DENSE SAND 
DISPLACED POSITION DENSE 

SAND 


GEOMETRIC SCALE 


DISPLACEMENT SCALE 


0 


Ft 
(3.05 mn) 
9.—Displacement Comparison 16-in. (406 Pipe Bend 


Based this analysis, design recommendations similar those shown 
Fig. developed for this study, should formulated which would allow tolerable 
bend movement and design conditions approaching fully restrained condition 
(Fig. 1). 

recognized that achievement high quality compaction not 
construction expedient especially the bends are located difficult access 
areas. Furthermore, achievement the desired compaction level closest 
the pipe which the most critical area very difficult. alternative 
compaction, selected sample concrete sand with 15% cement added 
and mixed with sufficient water produce pumpable slurry. The resulting 
slurry was left stand for least days triaxial test molds that were 
about in. (152.4 mm) high and with 3-in. (76.2-mm) diam allow the removal 
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excess water. Subsequently, these samples were tested manner similar 
the testing dense sand (Fig. 4). The resulting failure envelope indicated 
the cement stabilized sand stronger than the dense sand. our judgement 
cement stabilized sand will exhibit bonding cohesion particles more 
than compacted sand and its angle internal friction will also higher 
thus making stronger overall material. Consequently, the dense-sand 
can replaced sand stabilized with about 15% cement. This alternative 
should construction expedient since sand, cement, and water can combined 
conveniently located mixing plant and pumped in-place means common 
grout pumps and attenuated pipes; compaction equipment will required 
the site the bends. should pointed out here that the cement stabilized 


GROUND SURFACE 


-8 
Active) Passive 
Zone 


42 in. 8000 Ib per fr. 6ft. 3.5ft. 1.75ft. -Oft. Dense Sand 
(1067 mm) (116,721 m per m) (1-83m) (1.07m) (0.53m) 4 


36 in. 4750 Ib per ft. 4.Sft. 3.0ft. 1.50ft. -Oft. Dense Sand 
( 914 mm) (69,303 n per m) (1.37m) (0.92m) (0.46m) .82m) 


16 in. 2300 Ib per ft 4.Sft. 3.0ft. 1.0ft. Sot Dense Sand 
( 406 mm) (33,557 nm per m) (1.37m) (0.92m) (0.31m) (0.82m) 


16 in. 3700 Ib per ft. 10.0ft. 3.0ft. 1.0ft. 2.7ft. Medium Dense 
( 406 mm) (53,984 nm per m) (3.05 ) (0.92m) (0.31m) -82m) Sand 


FIG. 10.—Design Recommendations 


sand will set into fairly hard mass and, the future pipe bends require 
uncovering, some breaking and chipping the cement stabilized backfill will 
required. For given project and conditions, various sand-cement ratios 
should experimented with determine the most cost effective mixture. 


When designing buried pipelines, the pipe-soil interaction force, required 
restrain bend segment should computed and compared with the available 
soil resistance calculated conventional methods. larger than available 
resistance, use finite-element procedure and nonlinear soil properties will 
result economical design. 

Use finite-element technique and the nonlinear soil properties permit 
delineation the zone requiring select backfill which smaller than suggested 
conventional theories. This can result design and construction cost savings. 

For given pipe-soil interaction force, increase the depth soil 
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cover the most practical method achieving the higher backfill resistance. 
compacted sand backfill can replaced pumpable cement-sand 
slurry for construction expediency without any loss the backfill strength. 
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The following symbols are used this paper: 


intercept Hyperbolic stress-strain curve; 
slope Hyperbolic stress-strain curve; 
modulus elasticity; 

resultant end force pipe bend; 
soil cover; 

modulus number; 

modulus exponent; 

internal pressure pipeline fluid; 
atmospheric pressure; 

pipe-soil interaction force; 
horizontal axis p-g diagram; 
vertical axis diagram; 
bend radius; 

pipe radius; 

pipe wall thickness; 

coefficient thermal expansion; 
temperature differential; 

strain; 

Poisson’s ratio; 

hoop stress; 

longitudinal stress; 
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total major principal stress; 

total minor principal stress; 
effective major principal stress; 
effective minor principal stress; and 
maximum shear stress. 
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SUPERPLASTIC CONE-SHAPED SANDWICH 


(Reviewed the Aerospace Division) 


The applications sandwich construction are increasing steadily. addition 
aerospace vehicle structures, sandwich panels are also widely used housing, 
transportation, and packaging industries take advantage the high strength-to- 
density and high stiffness-to-density ratios. Traditionally, sandwich cores appear 
hexagonal cell shape and with vertical walls, which are often referred 
honeycomb. The use aluminum foil manufacture such honeycomb 
practically dominates all important structural applications sandwich panels. 
There are several disadvantages associated with the manufacturing and application 
such honeycomb cores: (1) involves expensive and complicated tooling 
facilities; (2) the core bonded facing sheets only the edge lines 
the core; and (3) the same geometrical configuration used for all applications 
and possibility provided for optimization studies. 

new method using superplastic alloy sheets and gas-forming techniques 
developed the Georgia Institute Technology (1) can eliminate all these 
drawbacks. this new process, different core configurations can obtained 
heating the superpiastic sheet critical temperature, applying gas pressure, 
vacuum pressure, both, and allowing the sheet collapse gently onto 
die having the desired geometry. Thus, through change modification 
the die, core different projection shape and positional symmetry becomes 
possible, which can studied experimentally and analytically for optimization 
purposes depending the load environment and application. using truncated 
projections, the top surfaces are available for bonding instead edges only. 
Shapes like hollow truncated cones, hollow truncated square hexagonal 

the October 22-26, 1979, ASCE Annual Convention and Exposition, held 
Atlanta, Ga. (Preprint 3694). 

School Engrg. Sci. and Mechanics, Georgia Inst. Tech., Atlanta, Ga. 
30332. 

Note.—Discussion open until June 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication November 
1979. This paper part the Transportation Engineering Journal ASCE, Proceedings 
the American Society Civil Engineers, ASCE, Vol. 107, No. January, 
1981. ISSN 0569-7891 $01.00. 
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pyramids, etc., arrayed square, hexagonal, etc., symmetry patterns, are just 
some examples from the possible configurations. Optimization studies undertaken 
include, e.g., the maximization the shear modulus the core (2,3,4,5), and 
the analysis configuration for least weight (6,7). this paper, procedures 
for the determination upper and lower bounds shear modulus for hollow, 


truncated, cone-shaped core (see Fig. with simple, square array symmetry, 
are reported. 


Properties 


The most important characteristic superplastic materials their exceptional 
stability uniaxial tensile deformation. Extremely large elongation 4,850% 
has been reported recently the literature for superplastic alloy Pb-62% 
eutectic without fracture necking, Langdon and his associates (8). This 
with most metals and metallic alloys where necking breaking usually occurs 
the range between and 50% elongation tension test. This extremely 
large elongation can take place low load and slow strain rate with temperature 
about half the absolute melting point without appreciable strain hardening. 


FIG. 1.—Sandwich Core Consisting Truncated, Hollow, Cone-Shaped Projections 


The essential microstructural condition extremely fine grain size, which 
under the correct conditions temperature and strain rate deforms largely 
grain boundary sliding and not the more usual slip and grain deformation 
processes. After deforming superplastically, relatively simple anneal restores 
the normal grain size and the normal metallic properties. 

The unusual superplastic property described here can provide simple and 
relatively inexpensive method producing complex shapes through various 
forming processes. has been reported that 30% savings weight and 50% 
cost are possible when superplastic forming combined with the diffusion 
bonding technique various aerospace structures (9). 


For efficient sandwich structures, the core must light weight and strong 
carrying the loadings, particularly the transverse shear. Therefore, the 
determination the shear modulus the core one the most important 
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tasks. order simplify the analysis, the following assumptions are made: 
(1) The complete core structure made identical projections, placed 
with square symmetry where the number such projections given panel 
large enough that any local, boundary effect, both, can ignored; 
(2) the effect the small radius the corners core projection due 
manufacturing negligible; (3) the wall thicknesses different locations 
the projections are the same, and are small comparison with other dimensions 
such the width, height, depth; (4) buckling warping the inclined 
walls the core projections excluded from consideration; and (5) the core 
assumed behave linearly elastic fashion. 

Kelsey, Gellatly, and Clark (10) made theoretical and experimental study 
the shear modulus foil honeycomb cores. Simple expressions for upper 
and lower limits were obtained application the unit displacement and 
unit load methods (11) conjunction with simplifying assumptions the 
strain and stress systems respectively, the core. Chang and Ebcioglu (12) 
studied similar problem. They adopted approach which corresponds 
the unit displacement method from equilibrium consideration. can proved 
that the results obtained Ref. correspond the upper limit solution presented 
Ref. 10. 

Using energy theorems (11), let the direction the applied force denoted 
subscript and the direction normal (in the horizontal xy-plane) 
subscript The shear stresses along directions and can expressed for 
the shear strains 


which the flexibility matrix. Conversely, one has 


which [k] the stiffnes matrix; and should noted that 
all practical tests the effective shear properties cores are limited the 
determination what effectively single direct flexibility i.e. 


Goce 

Upper Bounds.—Applying the principle virtual work the true state, 
and the virtual state, one obtains 


| {o,} "{€,} dv 


which k,, the required force applied produce unit displacement 
while all other displacements are zero. The variable represents the 
actual stresses due unit displacement with all others zero; the variable 
represents the virtual strains due unit displacement and the superscript 
denotes the transpose. should noted here that the 
used here are the unit overall shear strains and the core whole. 
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For the hollow, truncated cone-shaped core with square symmetry (see Fig. 
2), the stiffness can expressed 


which d,, and are parameters describing the dimensions the 
core; the shear modulus the core material bulk; and y,, the shear 
strain the wall due unit strain the core whole. 

Taking into account only the shear stresses the plane the inclined wall, 


FIG. 2.—Unit Projection Truncated, Hollow, Cone-Shaped Core 


the shear strain component inclined wall due unit overall shear strain, 


Similarly, the shear strain component the inclined wall due unit overall 
shear strain, the core the direction 


Substitution Eqs. and into Eq. yields 


Then from Eq. the effective shear modulus the core either xz- 
yz-plane the unit displacement method, can expressed 


d 


d = 


Lower Bounds.—The flexibility coefficients can expressed similar 
way 


which analogous Eq. The lower bound the effective shear modulus 


SANDWICH CORES 


can obtained using the unit load method. The ‘‘unit used here 
are the unit shear stresses and acting the core projection. The resulting 
shear flows developed the curved wall the projections due and 
computed from the conditions equilibrium within the core. Assuming 
that the shear stress constant for points the same elevation over the curved 
wall, can shown that f,, for core this particular shape can written 
the form 


sin 


which cos (see Fig. 3); and the shear flow the 
curved wall and elevation, sin from the top area the core 


FIG. 3.—Side-View and Sections Unit Projection 


due unit shear stress, acting the core. These shear flows can 
expressed 


sin 
Ecos 
Substitution Eq. into Eq. yields 


and 


| 
SECTION A-A 
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FIG. 4.—Nondimensional Effective Shear Modulus Versus a/b 0.5 
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FIG. 5.—Nondimensional Effective Shear Modulus Versus a/b 1.0 
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FIG. 6.—Nondimensional Effective Shear Modulus Versus a/b 1.5 
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FIG. 7.—Nondimensional Effective Shear Modulus Versus a/b 2.0 
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The effective shear modulus for the same core xz- yz-plane can now 
calculated the unit load method 


Anatysis 


The variables and given and 13, offer upper and 
lower bound respectively, for the effective shear modulus the core 
xz- yz-plane. They are programmed and plotted versus the angle shown 
Figs. for representative cases which 1.0, 1.5 and 2.0, 
respectively. The discrepancy between the upper and lower bound curves 
very small. can therefore concluded that the energy principles used 
this study are suitable for determining the shear modulus the sandwich core. 
The ratio G,/G, always less than unity, except for the case which 
This result follows naturally, since the methods used lead 
underestimate, and overestimate the true modulus. the special 
case, 7/2, the equality the two results indicates that the statically- 
equivalent stress system also kinematically compatible, that the two methods 
yield identical result and the analysis exact. From the representative curves 
shown, can observed that when the inclination angle equal 
the shear modulus always reaches its maximum value. For example, let 
value for approximately equal 1.40. Assuming Titanium 
4%V superplastic alloy used, with bulk shear modulus GPa and bulk 
density 4,500 kg/m’, then the effective shear modulus the core, 
turns out 483 MPa, which corresponds core density kg/m’. 
interesting observe that honeycomb core aluminum foil with 
density ft) can only offer shear modulus 228 MPa 
(33,000 psi) from Ref. 13. other words, one can gain factor two for 
the shear modulus the core the method presented this paper. 

The methodology developed here seems quite effective and useful for 
providing good engineering solution the problem. The new superplastically- 
formed core also offers competent approach for those who seek more efficient 
metallic structural design. combined weight and cost savings can therefore 
achieved, which important for the case which the structural 
weight critical factor. 
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The following symbols are used this paper: 


one haif the shortest distance between two neighboring unit projec- 
tions; 
the upper circle the truncated cone; 
length inclined wall; 
flexibility matrix; 
shear modulus bulk material; 
effective equivalent shear modulus core; 
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shear moduli the core, unit-displacement method and unit-load 
method, respectively; 

stiffness matrix; 

base length unit projection; 

thickness the wall; 

inclined angle the wall; 

shear strain; 

virtual strains due unit displacement and 

actual stresses due unit displacement 


S = 
B = 
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IMPACT BEHAVIOR BARRIERS 
NONLEVEL 


(Reviewed the Highway Division) 


AND 


Impact behavior roadside median barrier dependent number 
factors, including size and spacing posts, size and mounting height 
rail beam, offset beam from posts, soil conditions, and roadside conditions 
between the edge the traveled way and the barrier. Little known about 
the effects this last factor although may have the greatest influence 
performance. general, barriers have been designed and tested for flat terrain 
conditions even though roadside and median barriers are commonly placed 
side slopes, depressed medians, ditch bottoms, etc. 

the inception this study, trips were made several states survey 
current barrier placement practices and solicit input from various state 
transportation personnel. was found that there are four basic conditions for 
which roadside and median barriers are typically placed nonlevel terrain. 
These are shown the photographs Figs. 1-4. First, barriers used shield 
bridge piers, overhead sign bridge supports, other rigid objects depressed 
medians side slopes are often piaced near the object the barrier 
design permits. many cases this places the barrier the side slope, 
shown Fig. Secondly, barriers used shield bridge abutments other 
rigid objects near the shoulder are often flared away from the shoulder and 

“Presented the October 22-26, 1979, ASCE Annual Convention and Exposition, held 
Atlanta, Ga. 
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FIG. Barriers 


FIG. 1.—Barriers Shielding Rigid Objects 
| 
he, 
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terminated. consequence portion the barrier placed the side 
slope, shown Fig. Thirdly, roadside barriers are sometimes placed 
sections high-fill areas. Typically, the roadside composed 
shoulder, then relatively flat sloped embankment (usually 6:1 slope) 
which may extend (6.1 laterally, and finally, relatively steep 
embankment (usually 2:1 steeper). this case, the barrier placed 
the 6:1 slope shield the steeper embankment. This shown Fig. The 
last condition, which not common the other conditions mentioned, 
involves median barriers used prevent crossover head-on accidents that are 
placed stepped depressed medians. This shown Fig. 

gain insight behavior typical barriers nonlevel terrain, two roadside 
barriers were evaluated through full-scale crash test program. These were 
the system and the G4(1S) system. (Designations are consistent with those 


FIG. 3.—Barrier Barn Roof Section 


contained Ref. 4.) The system consists three steel cables mounted 
‘‘weak’’ posts, while the G4(1S) system consists standard steel W-beam 
mounted posts. Selection these systems was based four 
factors: (1) Both are operational (as per Ref. 1); (2) both are widely used; 
(3) both are used varying degrees nonlevel terrain conditions; and (4) 
the barrier system which may undergo considerable lateral 
movement impact, while the G4(1S) barrier system and 
will displace much less than the system for given set impact conditions. 

Selection 6:1 slope was based evaluation previous research 
involving computer simulation vehicle trajectory side slopes (see Appendix 
Ref. 4). The reference work showed that vehicle could possibly vault 
typical roadside barrier placed near the shoulder slopes steeper than approx 
10:1 8:1. was concluded that slope would the critical range 
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for lateral barrier offsets approx (3.66 from the edge the 
shoulder. 


FIG. 4.—Median Barriers 


This paper was condensed from report the Federal Highway Administration 
(FHWA) (6). 
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Test Site.—The test site, shown Fig. was designed represent common 
roadside section defined the road surface itself, shoulder, side slope, 
and back slope. The road surface consisted existing concrete surface 
(airport apron). 8-ft (2.4-m) shoulder area was cut 15:1 slope, and 
12-ft (3.7-m) side slope was cut 6:1 slope. (0.9-m-1.2-m) 
ditch bottom was added aid drainage and 4:1 back slope returned the 
configuration the original ground level. The existing horizontal concrete surface 
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FIG. 5.—Test Site Geometry 
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has downgrade slope parallel the rail and the test slopes were constructed 
parallel this drainage pattern. 

The length the installed roadside barrier was 200 (61.0 each case, 
excluding the rigid anchors both ends. Two sets anchor foundations were 
installed, one locating the face the guardrail (1.8 from the edge 
the shoulder (installation and the other locating (3.7 from 
the shoulder (installation B). 


FIG. 7.—G4(1S) Barrier 6:1 Slope 


| 


Description Barrier Systems.—Reference should made Fig. for details 
the two barrier systems described the following. 

G4(1S) Blocked-Out W-Beam.—This system consisted standard 12-gage 
steel W-beam mounted 8.5-steel posts. The blockout was 8.5, 
14-in. (35.6-cm) long steel block which was bolted the steel post. The top 
the rail was in. (68.6 cm) above ground level the post. The posts 
were spaced 6.25 (1.91 center-to-center, and they were set in. (106.7 
cm) 18-in. (45.7-cm) diam drilled hole and then backfilled with well-graded 
base material. This roadside barrier system, which was used tests 
shown Fig. 

Cable Guardrail.—This system consisted three 3/4-in. (1.9-cm) diam 
pretensioned cables mounted 5.7-steel posts. The top cable was 
in. (76.2 cm) from ground level, and the cables were spaced in. (7.6 cm) 
apart. The cables were attached the steel posts 5/16-in. (0.79-cm) diam 
steel hook bolts. 1/4 in. in. in.- (0.64 20.3 61.0 cm- 


FIG. Cable Rail Installed 6:1 Slope 


steel bearing plate was welded the back flange each post. The posts 
were spaced (4.9 apart center-to-center, and they were driven 29.5 
in. (74.5 cm) into the soil. 

Cables were attached the downstream end anchor turnbuckles, and 
the upstream end each cable was attached spring-compensating device. 
The cables were pretensioned approx 1,000 (454 kg) before the test. This 
system, used tests and shown Fig. 

Evaluation Criteria.—Ref. states that three appraisal factors are 
considered for the test longitudinal barrier: (1) Structural adequacy; (2) 
impact severity; and (3) vehicle trajectory hazard. summary, criteria for 
each appraisal factor are: 

Structural Adequacy.—The test article shall redirect the vehicle; thus the 
vehicle shall not penetrate vault over the installation. 

Impact the test article functions redirecting the vehicle, 
the maximum vehicle acceleration (50 m/sec average) measured near the center 


= 
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mass should less than the following values: 


Maximum Vehicle Accelerations (g’s) 
Lateral Longitudinal Total Remarks 


Preferred 
Acceptable 


These rigid body accelerations apply impact tests 15° less. 

Vehicle Trajectory Hazard.—After impact, the vehicle trajectory and final 
stopping position shall intrude minimum distance into adjacent traffic lanes. 

Two tests are recommended (5) evaluate roadside barrier along its 
need’’ (excludes end treatments). strength test which the barrier 
impacted mph (96.5 km/h) (2,043-kg) vehicle encroaching 
25°, recommended. The primary purpose this test evaluate the 
structural adequacy the barrier. The second test involves mph (96.5 
km/h) impact (1,022-kg) vehicle encroaching 15°. This test 
designed primarily evaluate the impact severity the barrier. 

Vehicle damage was assessed after each test accordance with two national 
rating scales, namely Traffic Accident Data scale (7) and the Society Automotive 
Engineers scale (2). 


Test 


Shown Table summary the six crash tests and results obtained 
therefrom. 

Test 1.—This was structural adequacy test the G4(1S) system for 
6-ft (1.8-m) offset. Upon impact, the lower part the bumper just cleared 
the top the rail and the car vaulted over the rail with little redirection. 

Test 2.—Test was nonstandard test but was conducted for two primary 
reasons: (1) determine the barrier was structurally adequate for lower 
encroachment angle (15° versus 25°); and (2) provide data for validating 
and updating barrier computer simulation programs (to used later supplement 
the crash test program). The car was contained and smoothly redirected. However, 
the car would have crossed adjacent traffic lanes, thereby posing potential 
trajectory hazard. 

Test 3.—This was structural adequacy test the G4(1S) system for 
12-ft (3.7-m) offset. Upon impact, the vehicle began redirect. However, the 
combined action lateral, longitudinal, and vertical loads the rail and right 
front wheel snagging support posts resulted complete fracture the 
rail. The right front wheel and portions the wheel assembly were torn free 
the vehicle. The vehicle penetrated the rail and came rest approx 
(27.5 beyond the point impact just behind the barrier. noted 
that post and rail material properties were compliance with recommended 
standards. 

Test 4.—This was impact severity test the G4(1S) system for 12-ft 
(3.7-m) offset. The car was contained and smoothly redirected. However, the 
car would have crossed adjacent traffic lanes, thereby posing potential trajectory 
hazard. 
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Test 5.—This was structural adequacy test the system for 6-ft 
(1.8-m) offset. The car was contained and smoothly redirected and there was 
trajectory hazard. Maximum dynamic deflection the cables was 9.5 
(2.9 m). 

Test 6.—This was impact severity test the system for 6-ft (1.8-m) 
offset. The car was contained and smoothly redirected, decelerations were below 


TABLE 1.—Summary Crash Test 


Adherence 
Performance 


Specifications 


speed, 
Barrier Vehicle Impact 
offset, weight, miles angle, 
(2) (4) (5) (6) 


“All barriers tested were placed 6:1 side slope. 

from outer edge shoulder face barrier. 

Section V-A for analysis performance specifications. 

vaulted over barrier. 

evaluation. 

penetrated through fractured rail element. 

Note: N/A Not Applicable; 0.305 0.454 kg; and mph 1.609 
km/h. 


suggested limits, and there was trajectory hazard. Maximum dynamic deflection 
the cables was 4.2 (1.3 m). 


Conclusions drawn result the six tests reported herein are: 


The G4(1S) roadside barrier system does not satisfy structural adequacy 
requirements when placed 6:1 slope offsets through (3.7 m). 
other words, the barrier, when placed stated, will not contain and redirect 
4,500 (2,043 kg) automobile impacting mph (96.5 km/h) and 
encroachment angle 25°. 

The G4(1S) system, when placed 6:1 slope and 6-ft (1.8-m) offset, 
will contain and smoothly redirect 4,500 (2,043 kg) automobile impacting 
mph and encroachment angle 15°. Although not proven the 
test, the writers’ opinion that the G4(1S) system will satisfy impact severity 


Struc- pact 
Test tural se- 
num- ade- |veri-| 
ber quacy? ty? angle? 
(1) (7) (8) (9) 
No* 
6.0 4,500 63.30 14.75 Yes 
12.0 2,300 58.2 14.75 Yes 
6.0 4,500 59.6 24.75 Yes |N/A Yes 
6.0 2,250 58.4 17.25 Yes Yes Yes 
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requirements when placed 6:1 slope and 6-ft (1.8-m) offset. 

The G4(1S) system satisfies impact severity requirements when placed 
6:1 slope 12-ft (3.7-m) offset. other words, the barrier, when placed 
stated, will contain and smoothly redirect (1,022-kg) automobile 
with tolerable decelerations when impacting mph (96.5 km/h) and 
encroachment angle 15°. 

Post-impact vehicle trajectory was less than desirable following the G4(1S) 
tests which the vehicle was redirected, i.e., tests and Results 
these two tests could interpreted mean that vehicle trajectory hazard 
existed; i.e., after impact, trajectory the vehicle would pose hazard 
traffic adjacent lanes. 

The roadside barrier system, when placed 6:1 slope 6-ft 
(1.8-m) offset, satisfied all three performance specifications for roadside barrier. 
When compared the G4(1S) system, improved performance the system 
attributed the 30-in. (76.2-cm) mounting height the top cable [versus 
in. (68.6 cm) for the W-beam]. The cables remained essentially the same 
height following impact, while the W-beam the G4(1S) system rotated backward 
and downward, creating ramp for the vehicle. 


RECOMMENDATIONS 


Results limited tests reported herein form basis for the tentative recommen- 
dations which follow. 


Roadside median barriers utilizing the standard W-beam, mounted 
in. (68.6 cm) above ground, should not placed 6:1 steeper slopes 
for offsets (3.7 m). Offset the lateral distance from the shoulder’s 
edge the face the barrier. When placed within these boundaries the barrier 
cannot expected contain and redirect automobile leaving the shoulder 
mph (96.5 km/h) with encroachment angle equaling exceeding 25°. 
other words, the barrier will not meet current performance specifications 
regarding (5). Barrier performance for offsets greater 
than (3.7 unknown and therefore recommendation can made 
regarding placement beyond this distance. Trajectory analysis (4) indicates that 
errant automobile with the aforementioned encroachment conditions would 
strike the barrier below the normal level terrain impact height for 
offsets greater than (3.7 m). must noted that the G4(1S) system 
did contain and redirect automobile 6-ft (1.8-m) and 12-ft (3.7-m) offsets 
for 60-mph (96.5-km/h), 15° encroachment angle. Statistics have shown that 
approx 80% all errant vehicles leave the travel way angle 15° 
less. Barrier systems now place 6:1 slopes utilizing 27-in. (68.6-cm) 
W-beam can thus expected contain and redirect large majority errant 
vehicles. 

The roadside barrier system acceptable for placement side slopes 
6:1 flatter. However, care must exercised its placement insure 
adequate distance behind the barrier for displacement during impact. lateral 
displacement approx (3.1 can expected the barrier impacted 
full-size car traveling mph (96.5 km/h) encroaching 25°. Therefore, 
the barrier should placed (3.1 more laterally from rigid objects. 


Tests (3) have shown the system will contain and redirect automobile 
when placed approx in. (45.7 cm) from embankments with slopes steep 
2:1 even though the barrier may deflect (3.1 more. However, 
the absence test data other supporting information, recommended 
that the system placed (3.1 more lateraily from embankments 
drop-offs with slopes steeper than 2:1. 
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SIMPLIFIED PROCEDURES FOR ESTIMATING 
ATTRACTIONS 


Ashok Kumar' 
(Reviewed the Urban Transportation Division) 


Estimates trips produced from and attracted traffic analysis zone 
are essential conventional four-step (trip generation, trip distribution, mode 
split, and route split) urban transportation modeling system (UTMS). This paper 
addresses the problems associated with obtaining reliable estimates trips 
attracted traffic analysis zone. Usually, rich data sets are available 
estimate the trips produced from traffic zone; e.g., origin-destination (O-D) 
surveys conducted home, Special Fourth Count United States Bureau 
Census Data Tabulations (24), census tract level population and household 
estimates prepared the National Planning Data Corporation (12), aerial 
photographs, etc. However, basic data pertaining variables needed estimate 
trips attracted traffic zone are usually either lacking difficult obtain. 

Variables commonly used estimate trips attracted zone are: (1) Gross 
floor area; (2) employment; and (3) gross land area. Statistical techniques such 
regression analysis and cross classification analysis correlate these variables 
trip attraction are well-documented two publications the Federal Highway 
Administration (7,20) and paper Shuldiner (16). The major problem 
estimating lies not the use these statistical techniques 
but assembling the data for the base year model calibration. The year model 
validation undertaken part the continuing 3-C process and the forecast 
year for which design traffic estimated. 

Transportation Planner, Northeast Ohio Areawide Coordinating Agency, 1501 Euclid 
Ave., Cleveland, Ohio 44115. 

Note.—Discussion open until June 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication November 14, 
1979. This paper part the Transportation Engineering Journal ASCE, Proceedings 
the American Society Civil Engineers, ASCE, Vol. 107, No. January, 
1981. ISSN 0569-7891 /$01.00. 
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Problems related data assembly the base year, model validation year, 
are mostly due budgetary constraints staff and data processing. some 
states the employment data type Standard Industrial Classification Code 
(SIC) and street addresses can obtained from the state bureau employment 
services the bureau workmen’s compensation, both. There are commer- 
cially available data sources such the listings compiled Dun and Bradstreet, 
but the completeness these listings remains seen. 

Despite the advances made setting the Dual Independent Map Encoding 
(DIME) File (5) for matching (ADMATCH) the street addresses census tract 
blocks (CTB) and subsequent compression these blocks traffic zones, 
very few large metropolitan areas have been able use the employment listings 
aforementioned estimate employment traffic zone. manual matching 
these listings traffic zones proves impractical, the study area 
large. Compilation gross floor area data (obtainable from Sanborn Maps 
and the county auditor’s office) the traffic zone level for the entire cordon 
area also cumbersome and impractical. 

The only readily obtainable data are the gross land area, secured from sources 
such the county auditor’s office; zoning maps prepared cities, townships 
and villages; United States Geological Survey (USGS) aerial maps and LANDSAT 
photographs (6). These data can also measured and summarized traffic 
zones inexperienced staff. 

Problems related data assembly for the forecast year are even greater. 
Large scale urban development models like Projective Land Use Model (PLUM) 
(2) and EMPIRIC Activity Allocation Model (2) are available the transportation 
planning literature for forecasting future year zonal population and employment. 
But, these models are not suitable for forecasting future growth and change 
for those metropolitan areas which are experiencing losses population and 
employment the central city. this writer’s knowledge, satisfactory 
mathematical models have been developed which can reflect decline, shift, and 
growth population and employment the traffic zone level for such urban 
areas. 

The Northeast Ohio Areawide Coordinating Agency (NOACA), the metropoli- 
tan planning organization (MPO) for the greater Cleveland area, faced similar 
problems when charged with the responsibility updating the areawide compre- 
hensive, land-use and transportation plans mandated the Department 
Housing and Urban Development, Federal Highway Administration, Urban Mass 
Transportation Administration and Environmental Protection Agency. The 
comprehensive transportation and land-use plan developed the 1960’s 
part the 3-C process using mathematical growth allocation model called 
Direct Trip Allocation Model (10,19) was found unsatisfactory replicating 
the base year (1975) traffic flows and land-use changes from 1963 (the year 
when O-D surveys were undertaken) 1975. 

was decided adopt different approach forecast future land-use changes. 
Local planners working together with the agency staff favored using process 
for allocating the future growth based upon the suitability analysis land available 


for Full documentation this procedure can 
found Refs. 13, 15. 
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ESTIMATING TRIP ATTRACTIONS 


Briefly, the method follows: available land judged the basis 
such attributes (1) Development site limitations like drainage problems and 
excessive slopes; (2) infrastructure availability like sewer and water connections; 
(3) site attractiveness like locale and quality educational system; (4) environ- 
mental limitations development like wetlands, poor-air noise 
zones; (5) transportation accessibility expressed terms nearness limited 
access highway interchange, mass transit facilities, and airport; and 
(6) development proclivity. 


TABLE 1.—Suitability Attribute Weightings Assumed Developing Trend Extended 
Land Use Plan 


Land Use Allocation Category 
High Low 


resi- resi- resi- 
Attribute dential dential dential 
(1) (2) (3) (4) 
Flood plains 
Poor drainage 
Steep slope 
Sewer service 
Water service 
Local alternatives 
Community service 
Protected areas 
Air/noise pollution 
Access interstate high- 
way 
Access dense arterial 
network 
Access rail 
Access airport 
High growth potential 
Industrial growth potential 


Note: Desirable attribute for development. Undesirable attribute for develop- 
ment. Restricted. present, attribute having weighting implies that future 
development that land use category cannot occur, even other desirable attributes 
are present. Required. absent, attribute having weighting implies that 
future development that land use category cannot take place, even other desirable 
attributes are present. 


The degree which attribute influences development varies from one 
land use another. After extensive examination with local planners, weightings 
for each the attributes for each land-use type were established (see Table 
1). series overlays were used map the attribute data. The area within 
the attribute outlines was shaded in. The attributes were mapped the negative 
direction; i.e., all positive suitability attributes were shaded where they were 
absent and all negative suitability attributes were shaded where they were present. 


Heavy 
in- in- Com- 
dustry dustry mercial 
(5) (6) (7) 
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This provided all-negative attribute maps. The reason for this procedure that 
after overlaying attribute maps, windows development are identifiable. 
computerized version this process now available through the Ohio Department 
Natural Resources (14). 

The future land-use allocation process can summarized follows: based 
upon historical data measured land-use inventories and the corresponding 
population and employment for the inventory years, relationship between 
population, employment, and acres consumed residential, commercial, and 
industrial land uses was established. Projected regional population and employ- 
ment changes 5-yr increments the year 2000 were converted into land-use 
acreage requirements. Using windows development identified suitability 
analysis, land-use acreage requirements were allocated. Finally, zonal boundaries 
were overlayed summarize the land-use changes 5-yr increments. 


Recognizing that assembling employment data the zonal level difficult 
not impossible, was decided investigate the possibility developing 
acceptable statistical relationships between person-trip attractions and the amount 
gross land area available for residential, commercial, industrial, and open 
space uses. Employment categories included calculating commercial and 
industrial land-consumption are shown Table 

was further decided express the amount residential land use number 
households rather than acres residential land use because compilation 


household data usually poses problem. was also decided express 
the amount commercial and industrial land use for zones falling within very 
high land-use intensity areas such the central business district (CBD) terms 
amount gross floor area rather than gross land area, because mixed 
land uses are found the same building. 


Several studies trip-generation analysis (e.g., Refs. 20), have indicated 
that the number trips generated from parcel land depends upon: (1) 
The amount; (2) the location; and (3) the type land use. Researchers have 
attempted identify the most prudent way express these variables correlating 
trips generated from given parcel land: the amount land (gross land 
area, gross floor area, number employees); location (CBD, central city, 
residential areas, distance from the CBD, regional accessibility); and type 
land use (manufacturing, retail trade, services, finance, insurance, and real estate, 
wholesale trade, transportation and public utilities, institutional land use). 

Studies done Black (4) indicate that the best statistical relationships are 
obtained when non-residential trip generation based upon floor area for 
commercial land use, upon employment for manufacturing (industrial) use, and 
upon land area for public buildings, public open spaces and transportation. 
Nakkash and Grecco (11) showed that using dummy variable indicate the 
location zone the central noncentral areas gave significantly better 
results than using gravitational-type regional access measure. Kolifarth and 
Shuldiner (9) showed that manufacturing land-use trip generation rates varied 
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TABLE 2.—Detailed Land Use Categories Under General Com- 
mercial, Retail, Office, and Open Space Land Use Variables Used Regression Analysis 


Industrial 
land use 
(1) 
Ordnance acces- 
sories 
Lumber wood 
products, except 
furniture 
Furniture fix- 
tures (household 
office furni- 
ture) 
Stone, clay and 
glass products 
Primary metal in- 
dustries 
Fabricated metal 
products 
Machinery, except 
electrical 
Electrical machin- 
ery, equipment, 
supplies 
Professional, scien- 
tific controll- 
ing instruments 
Transportation 
equipment 
Miscellaneous in- 
dustries 
Mining and quarry- 
ing 
Food kindred 
products (meat, 
dairy, etc.) 
Tobacco products 
Textile mill prod- 
ucts 
Apparel other 
finished products 
made from fab- 
rics 
Paper allied 
products 
Petroleum refining 
related indus- 
tries 


Commercial 
land use 


(2) 


Personal services 

Business services 

Automobile repair 
service 

Miscellaneous 
repair 

Financial, in- 
surance real 
estate 

Medical health 
services 

Other profes- 
sional services 

Eating drinking 
establishments 

Gasoline service 
stations 

Indoor theaters, 
auditoriums, 
field houses 

Other indoor rec- 
reation centers 

Wholesalers with 
stocks ware- 
houses 

Storage ware- 
houses 

Yard storage 

Raiiroad terminals 

Railroad yards 

Loca! interur- 
ban mass trans- 
portation sta- 
tions 

Motor freight 
transpor- 

tation—termi- 
nals 

Water transporta- 
tion—terminals 

Air transporta- 
tion—terminals 

Pipeline transpor- 
tation substa- 


Retail 
land 


(3) 


Shopping cen- 
ters 

Super markets 

Other food 
stores 

Drug stores 

Department 
stores 

Other general 
merchandise 
stores 

Apparel ac- 
cessories 
stores 

Furniture, 
home fur- 
nishings 
appliances 

Motor vehicle 
vehicle 
accessories 


terials 
Other retail 
stores 


Office space” 
(4) 
Office build- 

ings—gen- 
eral uses 
Federal gov- 
ernment— 
offices 
buildings 
State govern- 
ment—of- 
fices 
buildings 
County gov- 
ernment— 
offices 
buildings 
City govern- 
ment—of- 
fices 
buildings 


Open space 
(5) 


Public parks 
zoological 
gardens 

Golf courses 

Outdoor drive- 
theaters 

Racetracks 

Football 
baseball sta- 
diums 

Fairgrounds 

Other outdoor 
amusement 
centers 

Other public 
open space 
not classi- 
fied pre- 
viously 

Agriculture 

Forests 

Fisheries 
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Rubber products 
(tires, tubes, 
etc.) 

Leather products 


Communication 
Utilities—substa- 
tions 


Miscellaneous— Educational in- 
nondurable 
pens, etc.) schools 

Hospitals 


Other medical 
health buildings 
Cultural centers 


“Retail land use classification used only estimating home based shop attractions 
(HBSHPA). For other trip purposes, such as, home based work (HBW), home based 
other (HBO), home based social recreation (HBSRC) and truck trip ends estimation, 
retail land use lumped with commercial land use. 


space classification used only estimating trip attractions the CBD. 
For other locations, office land use lumped with commercial land use classification. 


significantly when manufacturing land-use types were stratified detailed 
two-digit standard industrial classification (SIC) codes (17) and the amount 
land use was expressed gross land area. However, when trip rates were 
derived using the number employees working within each detailed manufactur- 
ing SIC codes, significant variation was found. 

Building upon the research findings aforementioned and keeping perspective 
the data limitations, was decided develop the relationships between person 
trip-attractions and gross land area stratified generalized land-use types and 
some locational variable. Stratification traffic analysis zones location 
essential developing trip-attraction rates land area. 

During the early phase this study, scatter plots between zonal home-based 
person work trip attractions and the amount industrial and commercial land 
use within that zone were developed. Scatter plots were also developed after 
stratifying the zonal data pseudo-employment density measure, where 
pseudo-employment density defined the ratio between zonal home-based 
person work trip attractions and the amount zonal non-residential acreage 
(industrial plus commercial acreage). Correlations were found much higher 
after stratifying the data pseudo-employment density. Details this analysis 
can seen Ref. 23. 

Trip generation rates function land area for different land-use types 
have also been compiled the Institute Transportation Engineers (21) and 
the Maricopa Association Governments Transportation Planning Office (22). 
The trip rates developed these organizations are useful specific trip-generator 
analysis, but cannot readily used for the conventional four-step modeling 
process. The reasons are: (1) Only vehicular trips ends are estimated; (2) 
differentiation made between home-based and nonhome-based trip orientation 
and trip purpose; and (3) trip rates are not stratified location type. will 
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shown later that trip rates vary substatially location within the study 
area. 

For this study, pseudo-employment density measures locational variable 
was selected after considering other measures such residential density (number 
households per residential acre), gravitational accessibility, and distance from 
the CBD. The Cleveland metropolitan area multi-nucleus region with 
well-defined CBD Cleveland, Lorain, Painesville, and Elyria. Therefore, 
simple locational variable such the distance from the CBD not applicable. 
Gravitational access measure was found give similar stratification patterns 
pseudo-employment density and residential density. However, gravational 
access measure much more difficult compute compared residential and 
employment densities and hard comprehend non-quantitative planning 
professionals and decision makers. The accessibility measure, therefore, was 
dropped from further consideration. Except for school person trip attraction, 
pseudo-employment density was found give much superior statistical correla- 
tions. For school person trip attractions, residential density was found 
better stratifying variable. 

Trip attraction rates were developed for seven different trip purposes, namely: 
(1) Home-based work attractions; (2) home-based other attractions (such 
personal business, medical, dental, eat meal); (3) home-based social recreation 
attractions; (4) home-based school attractions; (5) home-based shop attractions; 
(6) nonhome-based trip ends; and (7) truck trip ends. For large metropolitan 
areas essential stratify the travel many trip purposes because 
trip length frequency distributions and mean trip lengths associated with these 
purposes are found statistically different. Combining some these trip 
purposes can bias the friction factors calibrated for gravity models and subsequent 
trip distribution. Stratification travel trip purpose also aids the selection 
causal variables statistical analysis. 

Table shows the relative distribution trip ends trip purpose and land-use 
type the trip destination for the Chicago region (1). This table clearly shows, 
expected, that school trips are predominantly attracted institutional land 
use, while social recreational trips can attracted residential, open space, 
and commercial land use. 

Using pseudo-employment density the stratifying variable, was necessary 
determine the range values this variable; i.e., strata over which trip 
attraction rates were statistically different. dummy variable technique covered 
Suits (18) especially useful this type analysis. Suppose one wishes 
test the hypothesis that truck trip rates based upon industrial land acreage 
vary significantly for suburbs (pseudo-employment density 0.1-20 em- 
acre), central city (20.1-40.0 employees 
acre), and the CBD (40.1+ employees/nonresidential acre). Note that these 
values are for illustrative purpose only. test this hypothesis, the following 
regression model estimated: 


which INDLAND zonal industrial acres; dummy variable having 
value equal when employment density 0.1-20 
acre, otherwise; dummy variable having value equal when employment 
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density 40.1+ employees/nonresidential acre, otherwise; and 
regression coefficients estimated. 

After and are estimated along with their statistic (7), the hypothesis 
that trip rates vary between suburbs and central city accepted rejected 
based upon the values the statistic associated with and found 
statistically significant, rates are assumed different between the 
central city and the suburbs. Similarly, examining the significance 
decision can made whether trip rates vary between the central city and 
the CBD. Note that for the density range 20.0-40.0 employees nonresidential 


TABLE Destination Land Use and Purpose 1970 Weekday Person Trips 


Social/| sonal 
recrea- 
Land use” tional 
(1) (7) 
Residential 
Manufac- 
turing 
Commercial 
Public build- 
ings 
Public open 
space 
Transpor- 
tation, 
commu- 
nication 
utilities 
mining 
Total des- 
tinations 
within 
region 
are thousands. 
Land use categories used Chicago Area Transportation Study. 
Note: Source—a summary travel characteristics—Chicago Area Transportation Study 
and Northwestern Indiana Regional Planning Commission, pp. 42, 1977. 


acre, and assume value zero the equation shown the foregoing 
and the value coefficient represents the trip rate for this density range. 

observed earlier, was necessary determine the strata pseudo- 
employment density over which trip attraction rates were different. order 
determine these ranges, zonal observations were allocated large number 
strata. The lower and upper bound each stratum was determined 
counting the number ‘of least observations falling within each stratum. 
Using the dummy variable technique described the foregoing, each stratum 
was checked against the adjoining strata for statistical differences trip attraction 


ness 
~ 
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rates. statistical differences were found exist, adjoining strata were 
collapsed into single stratum. The process was repeated until all strata were 
found statistically different. 

The computer program UREGRE (25) issued Urban Mass Transportation 
Administration was used creating the dummy variables and estimating the 
regression coefficients. The results obtained for each trip purpose are summarized 
the following. The data set used for these runs was the 1963 O-D survey 


conducted the Ohio Department Transportation for the Cleveland 
Metropolitan Area (3). 


Work Person Trip (HBWA) 


Using the linear regression analysis technique with dummy variables, the 
following equation was derived for traffic zones falling outside the Cleveland 


TABLE 4.—Values Calibration Coefficients Used Eq. 


Pseudo-employment density range, 
Coefficient Coefficient employees per nonresidential acre 


(1) (2) (3) 


within parentheses indicate the values statistic associated with regression 
coefficients shown. 


1.06 0.84 
(2.2) 
2.92 2.17 1-2 
(24.9) (1.9) 
5.13 5.39 2-3 
(17.2) (6.3) 
7.96 9.79 
(21.7) (44.5) 
11.79 11.26 
(21.6) (9.2) 
17.61 13.17 6.5-10 
(34.9) (22.2) 
24.29 23.05 10-13 
(58.9) (38.2) 
29.75 25.72 13-16 
(40.4) (32.3) 
38.19 33.80 16-20 
(49.5) (52.2) 
45.64 45.98 20-25 
(110.1) 
55.45 52.04 25-30 
(54.5) (31.4) 
66.46 70.44 30-40 
(58.6) (60.1) 
93.26 109.88 40-54 
(59.7) (50.7) 
149.22 128.02 54-80 
(70.9) (17.4) 
221.47 175.46 80+ 
(59.3) (9.9) 
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CBD. Values calibration coefficients and for different pseudo-employment 
density ranges are shown Table 


number observations 894, value 2,903.1, 
error estimate (HBWA) 166.6 


which HBWA zonal home-based work person trip attractions; COMMLAND 
zonal gross commercial land area, acres; INDLAND zonal gross industrial 
land area, acres; Pseudo-employment density (zonal HBWA)/(2. nonresi- 
dential area, acres); and nonresidential area zonal gross commercial land 
area zonal gross industrial land area. 


Person Trip Attractions (HBOA) 


This trip purpose category includes trips made for personal business, eating 
meals, and medical/dental visits. The equation giving best overall statistical 
fit follows: 


0.91, number observations 902, value 725.8, 


error estimate (HBOA) 270.8 


which HBOA zonal home-based other person trip attractions; number 


TABLE 5.—Values Calibration Coefficients Used Eq. 


Pseudo- 
employment density 
range, employees per 
Coefficient Coefficient Coefficient nonresidential acre 


(2) (3) (4) 


“Numbers within parentheses indicate the values statistic associated with regression 
coefficients shown. 


(1) 
0.40 0.66 905.2 
(27.9)* (2.1) (14.9) 
0.40 2.11 905.2 1-3 
(27.9) (11.7) (14.9) 
0.40 7.12 905.2 3-6.5 
(27.9) (15.2) (14.9) 
0.40 9.86 905.2 6.5-16 
(27.9) (17.8) (14.9) 
0.40 18.24 1,052.6 16-30 
(27.9) (16.7) (17.2) 
0.40 22.64 1,052.6 30-40 
(27.9) (10.5) (17.2) 
0.40 27.34 1,554.8 40-54 
(27.9) (9.7) (20.6) 
0.40 41.05 1,554.8 54+ 
(27.9) (16.1) (20.6) 
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households the analysis zone; and dummy variable having value 
equal the analysis zone contains major hospital well defined 
suburban town center, both, which can satisfy eat meal personal business 
trip purpose, otherwise. 

Values calibration coefficients and for different pseudo-employment 
density ranges are shown Table The coefficients and were also estimated 
without incorporating the dummy variable Eq. The statistical fit obtained 
was much inferior 0.75) compared that obtained when included. 
This indicates that major traffic generators tend have higher trip attraction 
because their regional draw and greater awareness the part the 
trip maker. fact, equations for estimating home-based social recreation trip 
attractions, shopping trip attractions, school trip attractions, and truck trip ends 
were all significantly improved when dummy variables indicating the presence 
major trip attractor were included the regression analysis. 


Recreation Person Trip Attractions (HBSRCA) 


For this trip purpose, the following equation was developed: 
0.90, number observations 911, value 1,031.2, 
error estimate (HBSRCA) 339.3 


which HBSRCA zonal home-based social recreation person trip attractions; 


OPEN SPACE zonal gross open space area, acres; and dummy 
variable having value equal the analysis zone contains one more 
the following generators: public golf course having large draw area; picnic 
area, recreational area, beach etc. having regional draw; YMCA, high school 


TABLE 6.—Values Calibration Coefficients Used Eq. 


Pseudo- 
employ- 
ment den- 
sity range, 
employees 
per nonres- 
Coefficient Coefficient Coefficient Coefficient idential acre 


(2) (4) (5) 


within parentheses indicate the values statistic associated with regression 
coefficients shown. 


0.84 1.31 0.44 1,047.8 
(6.3) (3.7) (28.6) 
0.84 5.14 0.44 1,047.8 
(25.1) (10.5) (3.7) (28.6) 
0.84 10.94 0.44 1,047.8 6.5-10 
(25.1) (17.0) (3.7) (28.6) 
0.61 10.94 5.51 1,047.8 10+ 
(34.2) (17.0) (10.8) (28.6) 


JANUARY 1981 


stadium baseball diamond having spectator stands, movie theater complex; 
otherwise. Value calibration coefficients and for different 
pseudo-employment density ranges are shown Table 


Person Trip (HBSCLA) 


Since most school trips are destined institutional land use (see Table 3), 
which only small subset the generalized commercial land use category, 
school person trip attractions cannot estimated directly function 
the commercial land-use category. School enrollment totals, however, can 
easily compiled the zonal level using data provided the state board 
education and colleges and universities. Consequently, school person trip attrac- 
tion rates were developed function zonal school plus college enrollment. 
Rates were further stratified residential densities capture the differences 
mode arrival, since schools located dense residential areas can easily 


TABLE 7.—Values Calibration Coefficients Used Eq. 


Residential density range, 
Coefficient Coefficient households per residential acre 


(2) (3) 


“Numbers within parentheses indicate the values statistic associated with regression 
coefficients shown. 

applicable—for residential density less than equal households 
acre, variable associated with coefficient not used 


accessed walking. The following equation was developed: 
number observations 590, value 551.8, 


error estimate (HBSCLA) 474.3 


which HBSCLA zonal home-based person school trip attractions; 
ENROLLMENT zonal school enrollment for grades kindergarten through 
college; dummy variable having value equal the analysis zone 
contains high school college which commuter oriented, otherwise; 
and residential density number households/gross residential land area, 
acres. 

Values calibration coefficients and for different residential density 
ranges are shown Table should noted that for residential density 


(1) 
1.41 
(36.8)* 
0.76 2-6 
(24.6) 
0.60 317.7 6-15 
(21.4) (2.6) 
0.29 798.1 15+ 
(7.7) (6.4) 
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values less than equal six households/acre, the dummy variable 
not used Eq. 

Several attempts were made improve the by: (1) Aggregating zonal 
data districts eliminate any bias caused due poor trip-file coding; (2) 
disaggregating the enrollment data different grades; i.e., kindergarten eighth, 
eighth high school, and college; and (3) separating models for school and 
college enrollment; but none these efforts proved useful improving the 
further. 


Home-based shop person trip attraction rates were first calibrated function 
the generalized commercial land-use category and stratified employment 


TABLE 8.—Values Calibration Coefficients Used Eq. 


Pseudo- 
employment density 
range, employees per 
Coefficient Coefficient Coefficient nonresidential acre 


(1) (2) (3) (4) 


11,990.0 
(41.4) 
11,990.0 
(41.4) 
11,990.0 
(41.4) 
11,990.0 
(41.4) 
11,990.0 
(41.4) 
11,990.0 
(41.4) 
within parentheses indicate the values statistic associated with regression 
coefficients shown. 


density, but considerably low value (0.26) was obtained. attempt was 
also made combine HBOA and HBSHPA into single trip purpose and 
estimate models similar the ones presented earlier for HBOA. This formulation 
also resulted low (0.40). 

became necessary base HBSHPA the retail land-use category rather 
than the generalized commercial land use. This puts additional burden 
base year data collection and design year forecasting land-use changes. 
Next work trips, shop trips usually account for the largest portion overall 
urban travel demand, this additional effort seems justified. Using retail land 
use the base variable, the following equation was developed: 


0.94, number observations 754, value 1,466.4, 
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error estimate (HBSHPA) 639.1 


which HBSHPA zonal home-based shop person trip attractions; RETAIL 
LAND zonal gross retail land area, acres; dummy variable having 
value equal the analysis zone contains subregional shopping center 
[characterized discount department store (s)/supermarket and group 
small specialty shops, drug store etc.], otherwise; and dummy variable 
having value equal the analysis zone contains regional shopping mall 
center (characterized two more nationally known department stores, 
many small specialty shops, and movie theaters), otherwise. Values calibration 
coefficients and for different pseudo-employment density ranges are 
shown Table 


Person Trip Enos (NHB) 


Since nonhome-based trips originating (O) from zone are found equal 
nonhome-based trips destined (D) that zone, model was calibrated 
estimate total nonhome-based person trip ends D). Estimated trip ends 
then can divided two and set equal nonhome-based trip origins and 
trip destinations. Another novel approach used estimate nonhome-based trip 
ends was correlate these trip ends home-based trip attractions rather than 
use number households and land area measurements land-use type. 
The following equation was developed for all zones falling outside the CBD: 


(25.98) (24.29) (17.58) 


933, value 4750.73, error estimate (NHB) 248.53 


which NHB nonhome-based person trip ends. 


with nonhome-based trips, truck trips originating from zone are found 
equal truck trips destined that zone. Therefore, models were calibrated 
directly estimate total truck trips ends. The following models were estimated: 


0.92, number observations 887, value 673.7, 
error estimate (truck trip ends) 202.9 


which dummy variable having value equal the zone contains 
truck terminal, industrial park, warehouse, established suburban town center, 
otherwise. Values calibration coefficients and for different 
pseudo-employment density ranges are shown Table 


Person Trip CBD 


was indicated earlier that for very high density zones such the CBD, 
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trip rates cannot based simply upon the gross land area because the 
mixed land-uses found the same building. Gross floor area was utilized 
the base variable. Models calibrated using the Cleveland CBD data are shown 
the following. Similar relationships have also been reported Wright (26). 


HBWA 3.522 OFFSPACE 5.876 COMMSPACE 0.449 
(3.8) (4.9) (0.3) 


INDSPACE; 0.82, number observations 35, value 46.7, 
error estimate (HBWA) 3,313.1 


HBOA 0.789 OFFSPACE 1.795 COMMSPACE; 
(2.5) (4.5) 


0.73, number observations 35, value 42.9, 

error estimate (HBOA) 1,145.7 

HBSRCA 1.058 COMMSPACE 1.382 RES SPACE; 
(6.9) (1.5) 


0.65, number observations 35, value 29.7, 


error estimate (HBSRCA) 599.4 


TABLE 9.—Values Calibration Coefficients Used Eq. 


employ- 
ment den- 
sity range, 
employees 
per nonres- 
Coefficient Coefficient Coefficient Coefficient idential acre 


(2) (3) (4) (5) 


“Numbers within parentheses indicate the values statistic associated with regression 
coefficients shown. 


Pseudo- 
(1) 
0.32 0.67 1.14 827.6 
(5.5) (2.7) (25.1) 
0.32 3.97 1.85 827.6 3-6.5 
(28.5) (9.9) (3.6) (25.1) 
0.32 6.00 1.85 827.6 6.5-10 
(28.5) (10.2) (3.6) (25.1) 
0.32 9.96 3.35 827.6 10-13 
(28.5) (12.8) (4.3) (25.1) 
0.32 12.41 6.27 827.6 13-16 
(28.5) (15.7) (13.9) (25.1) 
0.32 14.77 6.27 827.6 16-25 
(28.5) (19.3) (13.9) (25.1) 
0.32 14.77 21.65 827.6 25-30 
(28.5) (19.3) (20.2) (25.1) 
0.32 23.44 21.65 827.6 30-54 
(28.5) (19.8) (20.2) (25.1) 
0.32 23.44 84.26 827.6 
(28.5) (19.8) (10.9) (25.1) 
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0.91, number observations 35, value 324.2, 

error estimate (HBSHPA) 1,554.3 

NHB 0.462 HBWA 0.422 HBOA 0.138 HBSHPA 0.736 
(6.6) (1.5) (2.9) (2.9 


HBSRCA 0.99, number observations 35, 

value 681.9, error estimate (NHB) 311.5 

Truck Trip Ends 2.055 INDSPACE 0.783 COMMSPACE 0.204 
(5.5) (4.6) (1.6) 


OFFSPACE; 0.80, number observations 35, value 42.6, 
error estimate (truck trip ends) 470.4 


which OFFSPACE zonal gross office floorspace area, thousands 
square feet; COMMSPACE zonal gross commercial floorspace area, 
thousands square feet; RETAIL SPACE zonal gross retail floorspace area, 
thousands square feet; INDSPACE zonal gross industrial floorspace 
area, thousands square feet; and RES SPACE zonal gross residential 
floorspace area, thousands square feet. 

should noted that home-based school person trip attraction estimated 
using student enrollment and trip rates presented earlier. Notice that for the 
CBD models, office floorspace has been categorized separately rather than being 
lumped with the general commercial floor use category. 


The results reported here indicate that acceptable estimates person trip 
attractions can obtained using readily available data such the number 
households, the school enrollment, the area industrial, commercial, open 
space, and retail land. For zones falling within very high density area such 
the CBD, trip attraction should based upon floorspace measurements 
rather than upon land area. Categorization floorspace office, retail, 
commercial, and industrial use should adequate for estimating the person-trip 
attractions trip purpose. 

Stratification zones locational variable also shown essential 
because trip rates based upon gross land area vary significantly from low- 
employment high-employment density areas. Use pseudo-employment 
density was found adequate stratifying the zones. Introduction 
dummy variable reflect the presence major trip attractor relevant 
the trip purpose under consideration significantly improves the statistical correla- 
tion between trips attracted and the gross land area present within zone. 

For the model calibration year, pseudo-employment densities can easily 
calculated from the O-D survey and land-use inventory data. The year model 
validation has undertaken part the continuing planning process, 
pseudo-employment densities for zones falling within the mature portion 
urban area can assumed have the same values the base year 
(model calibration year). 
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For traffic zones belonging that portion the urban area which underwent 
new adjustment pseudo-employment density 
values can made two ways. One approach based upon the assumption 
that zone was partially developed (about 60-70%), the new development 
would follow the density that prevailed the base year. zone was 
underdeveloped, employment density can estimated averaging the prevailing 
densities neighboring zones. 

The second approach suitable the precise nature new development 
known. The analyst could either refer the comprehensive tabulation prepared 
Ide (8) which gives employment densities for establishments classified 
two-digit SIC codes, consult Ref. 23. Densities shown Ref. are based 
upon the pseudo-employment densities observed northeast Ohio. These sources 
can also aid the selection proper equations for assessing traffic impact 
associated with specific land development proposals. 

For long-range planning purposes, suggested that pseudo-employment 
densities used recursive manner, that is, base year densities used 
for estimating trip attractions for the first time increment (usually yr). Calculated 
work attractions then should used calculate the revised pseudo-employment 
densities, which turn should used for calculating trip attractions for the 
next time increment and on. 

Trip-estimating equations presented this paper can aggregated produce 
total person trip attraction. modal split and auto-occupancy these trips 
can estimated, vehicular trip generation estimates similar those reported 
the Institute Transportation Engineers (21) and the Maricopa Association 
Government (22) can obtained. should emphasized that certain traffic 
generators such airports, military bases, post offices, stadiums, large universi- 
ties, hospitals, museums, and libraries may require special adjustments 
independent trip-generation estimates. should also noted that trips calculated 


using the equations presented this paper account for travel occurring entirely 
within the cordon line. 
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The following symbols are used this paper: 


COMMLAND zonal gross commercial land area, acres; 


1 
19 


TE1 
COMMSPACE 


ENROLLMENT 


HBOA 
HBSCLA 
HBSHPA 
HBSRCA 


HBWA 

INDLAND 
INDSPACE 


NHB 
OFFSPACE 


OPEN SPACE 
RES SPACE 


RETAIL LAND 
RETAIL SPACE 
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zonal gross commercial floorspace area, thousands 
square feet; 

zonal school enrollment for grades kindergarten through 
college; 

zonal home-based other person trip attractions; 

zonal home-based school person trip attractions; 

zonal home-based shop person trip attractions; 

zonal home-based social recreation person trip attrac- 
tions; 

zonal home-based work person trip attractions; 

number households zone; 

dummy variable used regression analysis; 

dummy variable used regression analysis; 

zonal gross industrial land area, acres; 

zonal gross industrial floorspace area, thousands 
square feet; 

zonal nonhome-based person trip ends; 

zonal gross office floorspace area, thousands square 
feet; 

zonal gross open space area, acres; 

zonal gross residential floorspace area, thousands 
square feet; 

zonal gross retail land area, acres; 

zonal gross retail floorspace area, thousands square 


feet; 

calibration coefficient; 
calibration coefficient; 
calibration coefficient; and 
calibration coefficient. 
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The writer would like thank the discusser for his interest the paper 
and for offering comments regarding hoop stress analysis. 


Errata 


The following correction should made the original paper: 


CRACKING CONCRETE 


performance concrete pavements with cracking major concern 
and field with broad experimental and observational background, expansive 
and complex theoretical development and often gulf exists between the two. 


1979, Gerald Mouser (Proc. Paper 14691). 

Engr., Mountain Fuel Supply Co., 625 Connecticut Ave., P.O. Box 1129, 
Rock Springs, Wyo. 82901. 

1980, Lufti Raad and Jose Figueroa (Proc. Paper 15146). 

1980, Robert Deen, James Havens, Assaf Rahal, and Vernon 
Azevedo (Proc. Paper 15273). 

Soil Engr., Wisconsin Dept. Transportatiou, Div. Highways, 3502 Kinsman 
Blvd., Madison, Wisc. 53704. 
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The authors are complimented for another effort bridge this gulf. 

The writer has several comments that are essentially questions. Buckling and 
blowups are attributed expansion and implied these can prevented 
with expansion joints. this time there prevalent theory that expansion, 
growth, joint impaction other phenomena coupled with gravational forces 
cause creep pavement slabs down hill with forces sufficient move bridge 
abutments. attentuate this many agencies cut form expansion space 
back bridge ends. these forces exist and the expansion forces set forth 
the authors combine, should not almost all blowups occur sag points 
having bridges such those that fill over culverts rolling terrain? 
this the authors’ observation? 

The comments D-cracking are interesting most researchers have pre- 
viously attributed aggregate other mix properties with explanation 
lack edge deterioration midlengths the slab. This course lends 
credence the hypothesis fatigue loading the concrete temperature 
expansion-contraction presented the paper. But this hypothesis valid, 
D-cracking would occur all least most the joints for similar slab 
lengths built with similar joints. this true the pavements observed 
the authors? Naturally, one variation that could affect this different expansive 
properties the different pavements. this the reason that some pavements 
have D-cracking while others not, then the aggregate quality again probably 
becomes the major factor. some pavements D-cracked and others did not, 
were the thermal expansion wet-dry expansion properties the mix compo- 
nents, the cured concrete the various pavements determined? 


Errata 


The following corrections should made the original paper: 


Page 176, Table heading for Col. Should read FUEL USED”’ 

Page 181, paragraph line Should read and 10% diesel automobiles 
instead and 10% diesel automobiles, 


“July, 1980, Antoine Hobeika and Frank Tamberrino, Transportation and Develop- 
ment Around the Pacific, pp. 172-183. 
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